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ABSTRACT 

Preliminary  Turbine  Design  for  Thermal  Cycle  Feasibility 

David  A.  Potter 
Submitted  to  the  Department  of  Ocean  Engineering  and  the 
Department  of  Mechanical  Engineering  on  10  May  1974  in  partial 
fulfillment  of  the  requirements  for  the  degree  of  Master  of 
Science  in  Naval  Architecture  and  Marine  Engineering  and 
Master  of  Science  in  Mechanical  Engineering. 

Mathematical  models  of  an  axial  flow  turbine  and  a  radial 
inflow  turbine  were  developed.  Utilizing  these  models,  a  com- 
puter program  was  written  such  that  given  a  thermal  cycle  and 
a  minimum  of  input  variables,  the  necessary  turbine  design  cal- 
culations are  performed.  An  extensive  output  is  provided  such 
that  the  designer  may  utilize  this  output  in  order  to  evaluate 
the  practical  feasibility  of  the  given  thermal  cycle. 

The  Zener  sea  cycle  is  used  as  a  test  case  with  axial  and 
radial  designs  examined  for  feasibility. 

The  computer  program  is  written  in  the  Fortran  IV  language 
specifically  to  be  run  on  the  Interdata  70  computer. 
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NOTATION 

Blade  cross  sectional  hub  area 

Flow  area 

Annulus 

Blade  chord  length 

Specific  heat  at  constant 
pressure 

Specific  heat  at  constant 
volume 

Hydraulic  diameter 

Gravitational  constant 

Enthalpy 

Moment  of  Inertia  in 
tangential  direction 

Moment  of  Inertia  in 
Axial  direction 

Conversion  factor 

Windage  loss 

Stator  loss 

Rotor  loss 

Clearance  loss 

Mass  flow  rate 

Mach  Number 

Bending  moment  in  tangential 
direction 

Bending  moment  in  axial 
direc  tion 


UNITS 

ft2 

ft2 

ft2 

ft2 
BTU/°F 

BTU/°F 

ft 
32.174  ft  lbm/lbf  sec' 
BTU/lbm 
ft4 

ft4 

778  ft-lbf/BTU 
BTU/lbm 
BTU/lbm 
BTU/lbm 
BTU/lbm 
lbm/sec 
.  Non-Dim 
ft-lbf 

ft-lbf 


n  Number  of  axial  stages  Non-Dim 

N  Specific  speed  Non-Dim 

0  Throat  opening  feet 

P  Pressure  lbf/in 

0,  Volumetric  flow  rate  ft  /sec 

R  Universal  gas  constant  BTU/lbrn   R 

"R,  Reaction  Non-Dim 

r  Radius  feet 

R  Reynold's  Number  Non-Dim 

RPM  Shaft  rotative  speed  Rev/min 

S  Spacing  feet 

s  Entropy  BTU/lbm  °R 

T  Temperature  Deg  Rankine 

t  Blade  height  feet 

U  Rotor  tangential  speed  ft/sec 

V  True  fluid  speed  ft/sec 

W  Relative  fluid  speed  ft/sec 

x  Axial  distance  feet 

c^  True  fluid  flow  angle  from  Deg 
axial 

a  Relative  fluid  flow  angle  Deg 
from  axial 

vi  Ratio  of  C  /C  Non-Dim 
<J                P   v 


"l  Primary  loss  factor  Non-Dim 

?. 


Primary  loss  factor  corrected  Non-Dim 

for  aspect  ratio 
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?      Primary  loss  factor  corrected 

'"'  for  AH  k   H 
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Non-Dim 
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?  Slip  factor  Non-Dim 

yj  Efficiency  Non-Dim 

0  Tangential  angle  Deg 

*•  Stage  loading  coefficient  Non-Dim 

yU  Static  viscosity  lbm/hr  ft 

2 

*0  Kinematic  viscosity  ft  /hr 

A  Differential  function  Non-Dim 

tf~  Solidity  Non-Dim 

<S^  Centrifugal  stress  lbf/in 

YT  Constant  3.1415927 

O  Density  lbm/ft 

(y  Flow  coefficient  Non-Dim 

M^  Work  coefficient  Non-Dim 

6J  Angular  rotative  speed  Radians/sec 

SUBSCRIPTS 

0  Stagnation  conditions 

01  Machine  inlet  stagnation  condition 

02  Machine  outlet  stagnation  condition 
o2s  Isentropic  outlet  stagnation  condition 
i  Individual  flow  path 

t  Total  flow  paths 

b  Tangential  direction 

X  Axial  or  meridional  direction 
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m  Mean 

h  Hub 

t  Tip 

f  Fluid 

tt  Total  to  total 

ts  Total  to  static 

mat.  Material 

max  Maximum 

Axial 

1  Rotor  inlet 

2  Rotor  outlet 
Radial 

1  Stator  inlet 

2  Stator  outlet 

3  Rotor  inlet 

4  Rotor  outlet 
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INTRODUCTION 

The  thermal  cycle  is  the  basic  tool  in  use  today  to  con- 
vert the  energy  contained  in  a  source  or  fuel  into  useful 
power.   Cycles  using  geothermal  energy,  cycles  powered  by 
solar  energy,  and  cycles  powered  by  the  thermocline  existing 
in  the  oceans,  are  just  a  few  of  the  proposed  cycles  being 
studied  as  possible  sources  of  power  for  the  future.   In  the 
search  for  new  and  plentiful  energy  sources  in  the  years  to 
come,  new  and  different  thermal  cycles  will  be  proposed. 

The  evaluation  of  a  given  thermal  cycle  as  to  its  prac- 
tical feasibility  involves  a  detailed  study  in  which  not  only 
the  energy  source,  but  the  hardware  needed  to  convert  that 
energy  into  useful  power,  are  investigated.   In  such  a  study, 
the  engineer  must  concern  himself  with  three  main  areas:  con- 
verting the  energy  source  to  heat,  transferring  that  heat 
energy  to  a  working  medium,  and  then  utilizing  the  increased 
energy  of  the  working  medium  to  produce  useful  work  or  power. 
It  is  with  this  last  facet  of  the  problem  that  this  thesis  is 
concerned. 


13 

BACKGROUND 
1.1   Thermal  Cycles 

Thermal  cycles  come  in  all  sizes  and  shapes.   They  are 
used  not  only  to  convert  heat  energy  to  useful  work,  but  to 
transfer  heat  energy  from  one  location  to  another,  as  in  a 
refrigeration  unit.   The  difference  between  the  two  cycles, 
other  than  hardware,  lies  in  the  fact  that  in  a  work-orienta- 
ted thermal  cycle  energy  is  supplied  in  the  form  of  heat  or 
fuel  and  work  is  the  end  result;  while  in  a  heat  transfer 
or  refrigeration  cycle,  work  is  supplied  as  an  input  and  the 
end  result  is  the  transfer  of  heat  energy.   While  both  types 
of  cycles  operate  under  the  same  principles,  it  is  the  former 
type  which  is  of  great  concern  in  the  search  for  means  to 
supply  society's  growing  energy  demands,  and  which  will  be 
the  only  one  considered  from  hereon. 

Thermal  cycles  of  the  heat  to  work  type  can  be  catego- 
rized as  either  open  or  closed  cycles.   In  the  open  cycle, 
the  energy  of  the  source  (fuel)  is  transferred  to  a  working 
medium  in  the  form  of  heat,  whereupon  this  increased  energy 
of  the  working  medium  is  utilized  to  produce  useful  work. 
At  the  conclusion  of  the  open  cycle,  the  working  medium  is 
exhausted  from  the  machine.   In  a  closed  cycle  the  process 
is  essentially  the  same;  however,  the  working  medium  is  re- 
tained in  the  cycle  after  producing  work  and  is  reintroduced 
to  the  source  where  it  again  increases  its  energy  state  with 
the  addition  of  heat.   There  are,  of  course,  many  variations 
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and  combinations  of  the  above  cycles  using  many  different 
working  mediums  at  many  different  temperatures  and  pressures. 
1.2   Turbomachinery 

One  problem  common  to  all  thermal  cycles,  open  or  closed, 
is  that  of  finding  an  efficient  means  to  convert  the  energy 
of  the  working  medium  into  useful  work.   Without  considering 
unproven  schemes,  such  as  magnetohydrodynamics ,  there  are 
two  basic  classes  of  machines  in  use  today  with  which  this 
can  be  accomplished:   machines  which  extract  the  energy  from 
the  working  medium  as  a  change  in  its  potential  energy,  and 
those  machines  which  extract  the  energy  as  a  change  in  the 
kinetic  energy  of  the  working  medium.   Without  discussing 
the  relative  merits  of  either  type  of  machine,  it  is  suffi- 
cient to  state  that  the  second  class  of  machine,  the  turbo- 
machine,  is  a  viable  alternative  and  one  which  should  be 
considered  in  any  thermal  cycle  feasibility  study. 

Turbomachines ,  or  more  specifically,  turbines  used  in 
thermal  cycles  today,  can  be  divided  into  two  categories: 
those  in  which  the  working  fluid  flows  through  the  machine 
in  a  nearly  constant  axial  direction  parallel  to  the  center- 
line  of  the  machine;  and  those  in  which  the  working  fluid 
flows  in  a  direction  radial  to  the  centerline  of  the  machine, 
either  inward  or  outward.   The  choice  as  to  which  type  of 
machine  is  best  for  a  certain  set  of  conditions  is  a  diffi- 
cult one,  depending  on  a  large  number  of  variables  external 
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to  the  cycle  itself,  such  as  weight,  cost,  and  operating 
load.   This  will  be  discussed  more  fully  in  following 
sections. 
1.3   The  Design  Sequence 

in  determining  the  technical  feasibility  and  practical 
value  of  a  given  cycle,  at  the  preliminary  design  stage,  the 
engineer  must  follow  an  iterative  process.   One  logical 
sequence  would  be  to: 

1)  Assume  a  working  medium  and  a  set  of  states  for 
the  medium  (most  commonly  the  energy  extraction 
process  is  done  in  a  gaseous  state). 

2)  Assume  certain  hardware  restraints  and  a  flow 
rate  for  the  working  fluid. 

3)  Design  the  required  hardware  (if  not  available 
commercially) . 

4;   Evaluate  the  necessary  hardware  and  if  not 

satisfactory  return  to  step  1. 
This  process  is  repeated  several  times  until  the  hardware 
conforms  to  the  engineer's  definition  of  feasibility,  or  until 
some  prechosen  parameters  are  optimized.   The  problem  here 
lies  in  the  fact  that  the  calculations  necessary  to  design  the 
hardware,  although  relatively  straightforward,  are  numerous 
and  time-consuming.   For  a  given  set  of  cycle  conditions, 
there  are  often  many  options  open  to  the  designer,  all  of  which 
will  result  in  different  hardware  and  all  of  which  will  meet 
the  design  criteria.   The  choices  between  axial  or  radial 
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turbines,  and  between  turbines  or  reciprocating  machines, 
are  good  examples  of  this.   There  is  a  need,  therefore,  for 
the  designer  or  engineer  to  quickly  translate  his  cycle  re- 
quirements into  specific  hardware  sizes  and  characteristics 
in  order  to  logically  choose  between  the  options  available. 

1.4  Optimization 

Feasibility  and  optimization  are  two  words  whose 
definitions  can  be  different  for  every  different  cycle  con- 
sidered.  In  a  given  situation  it  may  be  necessary  to  max- 
imize efficiency,  while  in  another  situation,  the  design 
may  be  constrained  such  that  weight  or  volume  or  cost  must 
be  minimized.   As  for  turbines,  in  some  applications  an 
axial  flow  unit  might  be  chosen  for  its  high  efficiency, 
while  in  others  a  radial  machine  with  its  lower  fabrication 
cost  might  be  best.   Feasibility  also  has  different  meanings 
in  most  different  contexts.   In  one  instance  it  might  mean 
technically  possible,  while  in  another  it  might  mean  finan- 
cially profitable.   In  order  to  answer  these  questions,  the 
engineer  or  designer  must  first  define  feasibility  and  then 
decide  upon  which,  if  any,  parameters  are  to  be  optimized 
within  the  feasibility  constraints.   Having  done  this,  he 
may  then  evaluate  the  hardware  designs  and  make  a  logical 
choice  between  competing  designs. 

1 . 5  Purpose 

within  this  framework  then,  the  purpose  of  this  thesis 
is  to  bring  together  the  various  concepts  of  turbine  design; 
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to  develop  a  mathematical  model;  and  to  computerize  that 
model  in  such  a  way  that,  given  a  minimal  amount  of  in- 
formation about  the  cycle  and  the  machine  constraints,  a 
set  of  turbine  hardware  parameters  will  be  provided  that 
meet  the  requirements  of  the  preliminary  design  stage. 
These  may  then  be  used  in  the  design  sequence  mentioned 
previously  to  determine  feasibility.   No  attempt  will  be 
made  to  optimize  within  the  computer  program,  but  rather 
to  provide  information  from  which  a  suitable  optimization 
criteria  can  be  evaluated.   For  this  reason  no  attempt  will 
be  made  to  choose  between  an  axial  or  a  radial  turbine. 
Both  will  be  developed  simultaneously  and  information  will 
be  provided  on  both  types.   The  program  will  be  capable  of 
handling  multiple  inputs  for  the  same  cycle  in  one  computer 
run. 

It  is  hoped  that  by  utilizing  this  program  the  design 
sequence  may  be  shortened  considerably.   This  will  be  true 
especially  if  use  is  made  of  other  computer  programs  to 
calculate  heat  exchanger  parameters.   In  doing  so,  the 
designer  will  be  able  to  study  more  completely  the  possibi- 
lities of  a  given  cycle  at  the  preliminary  design  stuge, 
and  thus  develop  a  more  suitable  final  design. 
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DEVELOPMENT  OF  THE  MATHEMATICAL  MODEL 
2.1   General 

A  design  philosophy  is  a  necessary  ingredient  to  any 
worthwhile  design.   This  philosophy  may  not  always  be  ex- 
plicitly stated;  however,  it  is  always  present  in  a  good 
design.   It  is  a  means  whereby  the  designer  is  able  to  make 
logical  and,  of  most  importance,  consistent  decisions  during 
the  course  of  the  design  sequence.   This  is  especially  true 
at  the  preliminary  design  stage  where  decisions  are  made 
which  constrain  the  course  of  the  design  throughout. 

The  philosophy  followed  in  compounding  the  mathematical 
models  described  in  the  following  sections  and  in  the  actual 
programming  itself  can  be  stated  as  accurate  and  simple. 
Both  of  these  qualities  apply  to  any  work  done  at  a  prelimi- 
nary design  stage.   Accuracy  is  required  because  the  designs 
have  now  passed  beyond  the  back  of  the  envelope  stage;  and, 
although  detailed  design  is  not  required  yet,  the  results 
must  be  accurate  enough  so  that  detailed  design  can  begin  by 
using  them  as  a  basis.   Simplicity  is  a  must  because  at  the 
preliminary  stage,  neither  the  time  nor  the  money  are  usually 
available  to  carry  out  the  extensive  calculations  of  a  de- 
tailed design.   That  is  not  the  purpose  of  the  preliminary 
design.   One  last  idea  is  also  followed  in  the  handling  of 
certain  relationships;  that  in  the  absence  of  accurate  and 
simple  relationships  to  cover  a  certain  situation,  traditional 
and  reasonable  approximations  are  assumed.   With  the  lack  of 
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specific  data  in  certain  areas  this  borders  on  an  educated 
guess,  but  then  again  this  is  what  the  designer  is  paid  for. 

As  to  the  development  of  the  models  themselves,  they 
can  be  divided  into  three  main  areas:  representation  of  the 
working  medium  itself,  representation  of  the  axial  turbine, 
and  representation  of  the  radial  turbine.   They  will  be  ex- 
plained in  that  order. 

It  is  assumed  that  the  machines  described  by  the  models 
can  be  built  of  conventional  materials  and  that  the  axial 
blading  has  no  twist  and  is  set  in  a  purely  radial  direction. 
It  is  further  assumed  that  internally  cooled  blades  or  rotors 
are  not  utilized. 

All  computations  are  in  English  units  with  the  notation 
as  explained  on  pages  8  to  11.   Calculations  are  carried 
through  to  a  precision  consistent  with  the  accuracy  of  the 
Interdata  70  computer  (six  significant  figures). 
2.2   The  Working  Medium 
2.2.1   Equation  of  State 

Two  assumptions  concerning  the  working  medium  are  basic 
but  should  be  stated  for  clarity.   It  is  assumed  that  the 
working  medium  is  a  fluid  and  that  it  enters  the  machine  in 
a  gaseous  state;  this  precludes  the  design  of  hydraulic  tur- 
bines.  The  secoiid  assumption  is  that  the  inlet  and  exit 
states  of  the  fluid  are  known  and  that  the  designer  has  the 
properties  of  the  fluid  at  these  states  available  as  input. 

In  order  to  accurately  represent  the  fluid  within  the 
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machine,  an  equation  of  state  for  the  fluid  had  to  be  picked. 
The  perfect  gas  model  was  chosen  over  the  other  equations  of 
state  for  its  flexibility,  its  accuracy,  its  wide  acceptance, 
and  its  compatibility  with  the  standard  one  dimensional  com- 
pressible flow  relationships.   It  is  therefore  assumed  that 
in  the  region  under  consideration,  the  fluid  can  be  repre- 
sented by  the  perfect  gas  model  and  that  the  effects  of 
velocity  can  be  represented  by  the  standard  one  dimensional 
compressible  flow  relationships. 

To  make  this  assumption  more  realistic,  it  is  necessary 
to  calculate  the  different  constants  in  the  region  under  con- 
sideration rather  than  use  traditional  or  tabulated  values. 
One  method  to  achieve  this  is  to  use  data  for  an  isentropic 
expansion  in  the  region  and  to  calculate  the  constants  from 
this.   Since  the  designer  already  has  calculated  the  point 
in  order  to  determine  his  outlet  conditions,  it  was  decided 
to  use  an  isentropic  expansion  across  the  machine  for  the 
reference  points.   Thus  from  points  1  to  2s, 

C   =  (h  -  -h  . )/(T  _  -T  . ) 
p      o2s   ol  '   o2s   ol 

Using  the  Gibbs  equation  and  the  perfect  gas  law,  the 

following  commonly  used  relationship  for  an  isentropic 

expansion  may  be  derived: 

(P  .  /P  .)  =  (T  ,  /T  j  ^/UH-D 
v  o2s   ol       o2s'  ol 

where  J*=C  /(J 
P   v 
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Solving  this  equation  for  Jf1  : 

V    =  ln(P  -  /P  J/(ln(P  _  /P  .  )  -  ln(T  -  /T  .)) 
o2s   ol  '      02s'  ol         o2s   ol 

Thus   C   =  C  /  ja 

and  R  =  C   -  G  . 
P     v 

2.2.2  Flov  Rate 

In  describing  a  particular  turbomachine ,  the  mass  flow 
rate  of  the  machine  has  a  large  effect  on  its  design.   For  a 
particular  cycle,  the  mass  flow  rate  to  achieve  a  given 
output  can  be  easily  calculated  using  the  steady  flow  energy 
equation.   It  may  not  be  advantageous,  however,  to  achieve 
this  output  on  only  one  machine,  but  rather  to  utilize  several 
machines,  connected  in  parallel,  to  achieve  the  given  output. 
The  mass  flow  rate  through  a  given  machine,  with  no  allow- 
ance for  leakage,  will  be: 

ra   ,.  .  ,   ,  =  m,    ,  /  §   units  in  parallel 
individual     total  v 

2.2.3  Compressible  Flow  Relations 

The  two  parameters  necessitating  the  use  of  the  one 

dimensional  flow  relations  are  the  calculation  of  iMach  number 

and  the  calculation  of  flov  area.   Using  these  relationships, 

and  given  the  stagnation  or  total  temperature  and  the 

velocity  of  a  fluid  at  a  point,  the  Mach  number  and  flow 

area  may  be  easily  calculated.   From  the  definition  of  total 

temperature  and  Mach  number: 

T  =  T   -  V2/2g  JC 
o        feo   p 
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M   =   V/0fgokT)2 


In    its    common    form: 


m.(To)VP0Af    =   M(argo/K)V(l  +  (y-l/2)M2)(y+l)/2(^-l) 
Thus        Af    =    (milTo)VP0)    x    f(3\go,R,M) 


Af   =    (ra.(T    )2/F    )    x    (flow   function) 


where : 


Flow   function   =    (l  +  (r-l/2)M2  )  (^+1  )/2  &'1 )  /iM    UgQ/li)2 

2.3   Axial  Turbine 
2.3.1   General 

The  first  step  in  the  design  sequence  involves  deciding 
the  general  configuration  of  the  machine.   Is  the  hub  radius, 
the  mean  radius,  the  blade  tip  radius  or  some  combination  of 
these  to  be  maintained  constant  throughout  the  machine? 
Wilson  (l)  recommends  that  the  mean  radius  be  maintained 
constant,  while  Carmichael  (2)  recommends  hub  radius  be  kept 
constant  at  the  preliminary  design  stage.   Constant  mean 
radius  was  chosen  for  the  reason  that  the  effect  on  turbine 
performance  between  the  two  is  small,  and  calculations  at  a 
later  stage  will  be  simplified  with  a  constant  mean  radius. 
Following  this  assumption,  mean  conditions  are  calculated 
first.   As  the  velocity  triangles  and  other  blade  parameters 
vary  by  stage  throughout  the  machine  and  as  the  extremes  are 
represented  by  the  first  and  last  stages,  only  the  first  and 
last  stages  will  be  calculated  in  detail  with  the  other 
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stages  assumed  to  lie  in  between  these  values.   Finally, 
machine  parameters  which  depend  upon  the  previous  calcu- 
lation, such  as  total  efficiency  and  machine  length,  will 
be  made.   it  is  also  assumed  at  this  point  that  the  machine 
is  a  full  admission  turbine  operated  under  a  constant  load; 
both  assumptions  are  consistent  with  stationary  power  plant 
operation. 
2.3.2   Required  Output 

By  a  process  of  elimination,  the  list  of  quantities 
which  is  sufficient  to  describe  the  turbine  at  the  prelimi- 
nary design  stage  was  narrowed  to  the  following.   Many  other 
quantities  could  have  been  included;  however,  all  major  para- 
meters are  included,  and  those  omitted  may  simply  be  calcu- 
lated from  those  given. 


1)  First  and  Last  stage 
2) 
3) 
4) 
5) 
6) 

7)   Machine 
8) 
9) 
10) 


Tip  Diameters 

Blade  Height 

Aspect  Radios 

Number  of  Blades 

Efficiency  (  Total  to  total) 

Velocity  Triangles 

Critical  Mach  number 

Efficiency  (  Total  to  total) 

Critical  Stresses 

Length 
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In  describing  the  velocity  triangles  for  a  given  stage, 
three  locations  for  both  the  stator  and  rotor  entrance  tri- 
angles are  described:  hub,  mean  and  tip. 
2.3.3   Input 

In  listing  the  factors  which  are  generally  used  to  des- 
cribe a  turbine  stage,  it  was  found  that  they  could  be  cate- 
gorized in  three  groups:  those  which  describe  the  mass  flow 
rate;  those  which  describe  the  work  accomplished  per  pound 
of  fluid  per  stage;  and  those  which  describe  the  performance 
of  a  given  stage.   Each  facet  is  necessary  for  the  preliminary 
design . 

The  mass  flow  rate  for  the  machine,  and  thus  each  stage, 
(since  the  stages  are  assumed  to  be  in  series  with  no  bleed- 
offs)  can  be  described  by  the  total  mass  flow  rate  for  the 
cycle  and  the  number  of  turbines  in  parallel,  as  described 
in  section  2.2.2. 

The  work  per  pound  of  fluid  per  stage  can  be  simply  re- 
presented also.   Making  the  assumption  that  the  work  accom- 
plished per  stage  is  equal  for  all  stages  (2),  the  individual 
stage  work  can  be  represented  by  the  total  work  or  enthalpy 
drop  for  the  machine  and  the  number  of  stages. 

The  parameters  used  to  describe  stage  performance  are 
not  so  straightforward.   They  can  be  classed  as  either  pri- 
mary or  secondary,  depending  upon  the  strength  of  their  in- 
fluence upon  performance.   Primary  factors  are  generally 
associated  with  the  geometry  of  the  velocity  triangles,  and 
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if  the  deviation  is  assumed  negligible,  the  geometry  of  the 
blade  itself.   Secondary  factors  have  a  lesser  effect  on 
efficiency  and  are  related  to  the  likes  of  clearance,  aspect 
ratio,  and  other  mainly  mechanical  factors.   In  order  to  main- 
tain a  simple  machine  only  primary  factors  will  be  specified 
as  an  input,  with  a  logic  developed  for  the  determination  of 
reasonable  secondary  factors  based  on  the  given  primary 
factors . 

As  primary  factors,  Carmichael  (2)  recommends  reaction, 
R,  and  inlet  flow  angle,  c<  .   Craig  and  Cox  (3)  prefer  to 
use  the  stage  loading  parameter,  t,  also  called  the  work 
coefficient,  and  the  flow  coefficient,'?.  7\   is  another  widely 
used  parameter  (4).   Comparing  the  definitions  of  these 
parameters, 

*  =  x  "  (Vei  +  Ve2>/2U 

Tan  (<*x)  =  V0/Vx 
V  =AUVe/U2 

<j>  =  vx/u 

K    =   U/(V91  "  V02> 
it  can  be  seen  that  there  are  only  four  independent  variables. 
For  reasons  explained  in  a  later  section,  one. of  these  can  be 
eliminated;  thus  only  three  independent  variables  will  des- 
cribe the  shape  of  the  velocity  triangles. 
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The  three  variables  chosen  for  use  here  are~Rl,U  ,  and  o<s  . 
Reaction  was  chosen  for  its  compatibility  with  the  Euler  work 
equation  which  will  be  used  later  to  solve  for  a  mean  radius. 
Given  this  mean  radius  and  an  h.PM,  the  second  variable,  U, 
can  be  specified.   The  last  variable ,  oC1  ,  was  chosen  as  it 
would  have  to  be  calculated  anyway  during  the  course  of  the 
design. 

The  input  for  the  axial  turbine  can  be  summarized  as: 

1)  Mass  flow  rate 

2)  Number  of  units  in  parallel 

3)  Number  of  stages  per  unit 

4)  Shaft  revolutions  per  minute 

5)  Reaction  at  mean  conditions 

6)  oC,    at  mean  conditions 

7)  Cycle  conditions  of  section  2.2.1. 
2.3.4   Mean  Conditions 

The  starting  point  for  the  mean  conditions  is  with  the 
simultaneous  solution  of  the  Euler  work  equation  and  the 
reaction  equation.   It  is  assumed  that  the  fluid  enters  the 
stator  only  in  an  axial  direction,  with  no  tangential  compo- 
nent.  This  is  a  common  assumption  (4)  and  is  made  for  several 
reasons.   One  argument  is  that  the  turbine  works  by  extracting 
energy  from  the  fluid;  thus  any  energy  left  in  the  fluid  at 
rotor  exit  (stator  inlet)  in  the  form  of  vorticity  (tangential 
velocity)  is  unavailable  for  extraction,  and  thus  contributes 
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to  a  loss.   Using  this  assumption,  the  Euler  equation  applied 

to  mean  conditions  across  the  rotor  becomes 

Ah   =  60  r  V„,  /g  J 
o      m  91m  &o 


n 

where  Ah   =  h  , -h  „ 
o     ol   o2 

The  reaction  equation,  also  applied  at  mean  conditions  across 

the  rotor  is  simplified  to 

H     =  1  -  YQ1  /2wr  . 
m         81m     m 

The  simultaneous  solution  of  these  equations  for  r   and  V_-. 

^  m  01m 

with  W=TfRPM/30    gives 

r      =    (Ah    g  J450/n^HPM2n(l-U))2 
m  o&o  m 

V61m   =   rTrRPM(l-TU/l5. 
m  m 

Substituting  in  the  appropriate  constants, 

r   =  1068.247  (Ah  /n(l-R  ))5/RPM 
ra  o       m 

and 

V91m  =  223.733lVAho(l-U)/n)i 

In  order  to  prevent  the  solution  from  expanding  at  "R  =1,  H» 

r  r  6      m   '   m 

must  be  limited  to  less  than  0.99. 

The  rotor  speed  can  be  calculated  as 

U,   =  U-   =  U   =Yi-RPM  r  /30. 
lm     2m     m  ra 

Using  the  trignometric  relationships,  Figure  1,  and  the 

©C,  which  was  given  as  an  input,  the  rest  of  the  velocity 

triangle  for  the  rotor  inlet  can  be  calculated. 


Axial  Turbine  Velocity  Triangles 


Rotor  Inlet 
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Stator  Inlet 


Figure  1 

V,   =  VQ1  /Sin(<*.  ) 
lm     91m'    v  lm 

V  =  V    /Tan(c<  ) 
x     91m       lm 

/?lm   =  ArcTan(Tan(o<lm)  (1-1/(2-2^))) 

V  ' =  V  /Cos (A  ). 

lm     x       lm 

Making  the  assumption  that  the  axial  velocity  is  constant 

throughout  the  machine  as  was  stated  before,  the  tangential 

velocity  at  rotor  outlet  is  zero;  thus  the  velocity  triangle 

for  the  stator  inlet  (rotor  outlet)  can  be  calculated: 

'J3~      =  ArcTan(Tan(oC     )/2(l-R,)) 
~ 2m  lm  m 
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\<-        =  V_   /Cos  (/a,  )Tan(oC  ) 
2r     91m'    v  <>m     v  lm' 


oe        -  V     =  0 
2m     02m 


V  .  =  V  .  =  V 
xl     x2     x 


V„   =  V  . 
2m     x 

The  total  gas  path  deflection  for  the  rotor  is  the  sum 

of  the  inlet  and  outlet  angles  (/^   +  /S^.  );  while  for  the 

&       lm   ^2m  ' 

stator  it  is  just  c*   ,  as  the  inlet  angle  is  zero. 
2.3.5   Stage  Conditions 

In  solving  for  stage  conditions  it  can  be  assumed  that 
the  extreme  values  are  represented  by  the  first  and  last 
stages,  with  all  other  stages  falling  in  between.   If  the 
variance  is  assumed  linear,  the  dimensions  and  velocities 
can  be  calculated  for  any  stage  with  reasonable  accuracy. 
This  assumption  also  allows  the  calculation  of  machine 
efficiency  as  the  mean  of  first  and  last  stage  efficiencies. 
Based  on  these  assumptions,  the  following  calculations  can 
be  made  for  the  first  and  last  stages. 

Using  the  one  dimensional  compressible  flow  relationships 

described  in  section  2.2.3,  the  flow  area  can  be  calculated  for 

a  given  velocity  and  stagnation  temperature.   The  flow  area 

and  annulus  area  are  related  as  shown  in  Figure  2. 

Area     =  Area„,   /Cos(o4). 
ann        flow 

Using  the  definition  of  annulus  area  and  the  mean  radius, 
the  blade  height  may  be  found. 


Axial  Turbine  Flow  Area 
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Machine 
Axis 


Figure  2 


2    2 
Areaann  =^ (rt  "rh  ) 

=  YT(rt-rh)(rt+rh) 


=YTt2r 


m 

Blade  height  =  t  =  Area    /2  r  ft 

b  ann     m 

where  r   =  (r  +r  )/2 

m      I   h 

The  tip  and  hub  radii  are  found  by  respectively  adding  or 
subtracting  one  half  of  the  blade  height  to  the  mean  radius. 

The  next  step  involves  making  the  assumption  that  a 
streamline  of  fluid  has  no  appreciable  curvature  in  the  radial 
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direction  as  it  moves  through  the  machine.   This  is  not  com- 
pletely correct,  but  for  small  wall  curvatures  it  can  be  con- 
sidered accurate  (l).   This  assumption  means  that  a  segment 
of  fluid,  as  shown  in  Figure  3,  must  be  in  radial  equilibrium, 

KAUIAL  EQUILIBRIUM  FORCES 


p  +  d_o.  dr 
dr 


Axis  of  machine  into  page 


Figure  3 


Summuation  of  the  forces  acting  on  this  segment  gives  the 


equation  for  this  condition  of  radial  equilibrium, 

2 
0  '  eo" 


JL_  dp.  =  VQ  /g/ir 
/>f    dr 


Using  the  Gibbs  equation  and  the  definition  of  total 
enthalpy  for  an  adiabatic  condition,  the  radial  equilibrium 
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equation  can  be  rewritten  as 

dh  -  Tds  =  1    dV  2  +  1      d(r2V  2). 
o    x     2       9  ' 

dr      dr    2g   dr      2g  r      dr 

Making  the  assumption  that  the  enthalpy  drop  is  uniformly 

distributed  over  the  blade  radially,  and  that  entropy  is  not  a 

function  of  radius,  the  result  is 

dV  2  =  -1   d^r2V  ,'2). 
x     2  9 

dr       r      dr 

At  this  point  one  of  several  distributions  for  V   or  Vrt 
r  x      0 

may  be  chosen;  however,  having  made  the  assumption  previously 

that  axial  velocity  is  constant  throughout  the  machine,  the 

above  reduces  to 

r  V   =  constant, 
w 

This  constant  may  be  calculated  from  the  mean  conditions 
and,  using  the  above  relationship,  the  tangential  velocities 
at  hub  and  tip  can  be  calculated.   The  blade  speed  is  easily 
calculated  at  hub  and  tip  as  shown  previously;  and,  using 
trignometric  relationships,  the  complete  velocity  triangles 
for  rotor  and  stator  are  easily  calculated  in  the  manner 
illustrated  in  section  2.3.3. 

The  choice  of  aspect  ratio  (blade  height/blade  chord) 
must  now  be  made.   As  this  is  a  secondary  variable,  section 
2.3.3,  a  means  must  be  derived  for  a  logical  decision.   A 
study  of  data  presented  by  Forrester  (5)  shows  that  the 
secondary  losses  attributed  to  aspect  ratio  effects  decrease 
with  an  increasing  aspect  ratio;  thus  a  large  aspect  ratio 
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would  be  advantageous  from  an  efficiency  standpoint.   On  the 
other  hand,  Horlock  (6)  states  that  the  bending  stresses  on 
the  blade  may  place  a  limit,  on  the  blade  chord  length  and 
thus  limit  the  feasible  aspect  ratio.   Any  aspect  ratio 
chosen  should  therefore  be  as  large  as  possible  without  ex- 
ceeding a  predetermined  stress  limit.   Carmichael  recommends 
a  limit  of  one  ton  per  square  inch  which  was  used  in  the 
following  derivation. 

In  order  to  simplify  the  calculations  necessary,  the 
fluid  was  modeled  as  incompressible  and  the  blade  shape  was 
simplified  so  that  the  resulting  criteria  could  be  applied  to 
the  general  case.   The  maximum  stress  was  assumed  to  occur  in 
the  rotor  at  the  root  of  the  blade  and  located  on  the 


Blade  Root  Stresses 


Figure  4 
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upstream  side.   By  considering  the  change  in  tangential  momen- 
tum across  the  blade,  the  bending  moment  in  the  tangential 
direction  (Figure  4)  can  be.  shown  to  be 

M0=/VfV  r<r  -  rh)(V  dr 
rh    ~ 

where  z  =  number  of  blades  per  stage. 

Integrating  this  equation,  and  substituting  the  Euler 

work  equation  and  the  definition  of  blade  speed  gives 

M„  =  mah  J(r,  -r,  )260  S, 
6      o    t   h      1 


(r  2-r,  2)  (2K)2  RPM  r,  n 
t    h  h 


Simplifying  and  substituting  in  the  constants  yields 

M„  =  591.209  riiAh  (r,-r,  )  S1 /RPM  r   r,  n. 
W  o   t   h    1       m   n 

Since  Tan  (/3  )    =   Forces.,  ,.    ,  .   /  Forces   ,.    ,  . 

'm  9  direction         x  direction 

the  bending  moment  in  the  axial  direction  can  be  described  by 

M  =  M„  Tan  ( ,0   )  . 
x    9       m 

For  a  particular  design  let  Kl  be  a  constant  representing 
Kl  =  mail   (r,  -r,  )  591.209/RPM  r,  r   n 

O         h  X)  XI      ill 

then         MQ  =  Kl  Sx 

M   =K1  S.  Tan  {/3  )  . 
x        1       /m 

For  the  blade  shape  shown  in  Figure  4,  the  moments  of 
inertia  about  the  respective  neutral  axes  can  be  shown  to  be 

Ig  =  c  Sx3/69.429 

I-  =  Sj  c3/72  . 
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For    stress  at  the  point  A, 

e/1-  =  30.86/c  S  2 
x  1 

x/lg  =  36/SlC2. 

The  bending  stress  is  then 

o~h   =  Kl(Tan(^)36/c2  +  30.86/c  S  )  . 

Using  Zweifel's  criteria  to  relate  S  and  c  (4), 

0.8  =  2  Cos2(^2h)(Tan(^h)  +  TanU^))  S1/c 

substituting  in  the  bending  stress  equation  and  solving  for 
chord,  using  the  one  ton  per  square  inch  stress  limit  men- 
tioned previously, 

c  =  (K2  f)? 

where  K2  =  Kl  32 . 174/(2000) (144) 

f   =   36   Tan(^lm)    +   77.15    Cos2(/?2h)    x 

(Tan(/?lh)    +    Tan(/?2h)). 

This  approximation  for  chord  is  based  on  a  stress  limit 
caused  by  gas  bending  forces;  however,  to  account  for  the 
fact  that  other  design  conditions  may  impose  a  lower  limit 
than  this  criteria,  a  maximum  limit  to  aspect  ratio  of  10 
will  be  imposed. 

Zweifel's  criteria  can  now  be  used  to  determine  spacing 
for  the  rotor  from  this  calculated  chord.   Assuming  the  same 
chord  for  the  stator,  spacing  can  be  calculated  there  also. 
As  shown  in  Figure  5,  the  throat  opening  can  be  figured  for 
both  the  rotor  and  stator. 
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Throat  Opening 


Axis  of  machine 


Figure  5 


°1  =  Sl  C°S^lJ 

°2  =  S2  Cos^2m) 

Knowing  the  spacing,  the  number  of  blades  for  rotor  or 

stator  is  simply 

Z  =  2yrr,  /  S. 
h 

t 

Using  the  hydraulic  diameter  of  the  throat  opening, 
again  for  both  rotor  and  stator, 

D  =  2  0t/(0  +  t) 
the  Reynolds  number  for  both  the  rotor  and  stator  can  be 
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calculated  as 

He2  =  D2   Vm2/S> 

Itel  =  Dl  Vml^- 
2.3.6   Efficiency 

The  determination  of  efficiency  for  axial  turbines  is  a 
much-debated  subject.   Craig  &  Cox  (3),  Smith  (5),  and  Ainley 
and  Matheson  (5),  all  correlate  their  data  using  the  work 
and  flow  coefficients  mentioned  previously;  while  Soderberg 
(4)  correlates  using  gas  deflection  angle.   In  either  case 
the  primary  loss  is  related  to  blade  profile  loss.   Secondary 
losses  are  attributed  to  many  factors,  some  of  them  being: 

1)  Clearance  loss 

2)  Reynolds  number  effects 

3 )  Aspect  ratio 

4)  Guide  loss 

5)  Runner  loss 

6)  Guide  gland  loss 

7)  Lacing  wire  loss 

8)  balance  hole  loss 

9)  wetness  loss 

10)  Disc  windage  loss 

11)  Partial  admission  loss 

12)  Annulus  wall  loss 

13)  Cavity  loss. 

While  including  all  of  these  factors  would  give  a  most 
accurate  efficiency  prediction,  the  accuracy  required  in  a 
preliminary  design  doesn't  require  the  inclusion  of  them  all. 

Soderberg' s  correlation  was  chosen  in  that  it  provides 
accuracy  consistent  with  the  intent  of  the  program;  yet  it 
is  simple  enough  to  be  easily  programmed  on  a  computer.   This 
correlation  takes  into  account  profile  loss,  clearance  loss, 
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aspect  ratio,  and  a  Reynolds  number  effect.   These  losses  are 
equated  to  a  fraction  of  the  kinetic  energy  available  to  the 
turbine  and  a  total  to  total  efficiency  calculated. 

The  method  consists  of  finding  a  primary  loss  factor, 
based  on  gas  deflection  angle  for  both  stator  and  rotor. 
This  loss  factor  is  modified  for  aspect  ratio  and  Reynolds 
number  effects,  and  then  an  efficiency  calculated.   The 
efficiency  is  then  modified  by  a  factor  which  takes  into 
account  the  clearance  losses. 

?   =  f  (Gas  Deflection) 


h  =  ((1+?o)(-975+0.075  c/t))-l 
?  -  7   (^} 


5  «25 


The  above  are  calculated  separately  for  both  stator  and  rotor. 

The  efficiency  is  then: 

2    -,   D  „    2x-l 


V  1    =    (W2Stator    Vl  +  ?r>Roto 


r  V  ') 


2U<W 

^?+  +  (Stage)  =  V),  (Area  Annulus  -  Tip  Clearance  Area) 

Area  Annulus 

The  primary  loss  coefficient  is  one  which  can  be  pro- 
grammed easily  in  a  curve-fitting  subprogram.   The  tip  clear- 
ance area  was  calculated  based  on  a  tip  clearance  of  one  per- 
cent of  the  tip  radius.   This  gave  a  constant  factor  of  0.99. 
Assuming  a  linear  distribution  of  efficiencies  between  the 
first  and  last  stages  means  that  the  machine  efficiency  is 
the  average  efficiency. 
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l   Machine  =^First  Stage  +  >^Last  Stage  ''2' 

2.3.7  Miscellaneous  Machine  Parameters 

In  order  to  more  easily  picture  the  machine  size,  all 
pertinent  radii  are  converted  to  diameters. 

A  machine  unit  length  is  calculated  using  the  average 
stage  chord  and  the  number  of  stages. 

The  maximum  blade  centrifugal  stress  occurs  in  the  last 

stage  at  the  rotor  blade  root.   This  can  be  written  as 

follows: 

<5~  =  /mat.     (        co  r  A,  dr . 

c    : —      J  h 

K  A,      ' 

h 

If  the  blade  is  assumed  to  have  no  taper,  A,  is  constant. 

h 

Substituting  in  the  constants  and  assuming  the  density  of 

steel  is  approximately  0.29  lbm/in   gives 

<T  =  0.001186  RPM2  r   t. 
c  m 

2.3.8  Turbine  Limitations 

Certain  limitations  must  be  placed  on  the  mathematical 
model  previously  described,  so  that  it  won't  produce  turbine 
parameters  for  a  machine  that  is  either  physically  impossible 
or  not  consistent  with  the  assumptions  of  the  model.   These 
limitations  will  be  in  the  form  of  numerical  checks  in  the 
program  and  will  cause  the  calculations  to  be  terminated  if 
they  are  not  satisfied. 

The  first  check  is  for  excessive  gas  deflection.   While 
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this  in  itself  is  not  an  impossible  situation,  only  leading  to 
a  lower  efficiency,  the  data  available  for  the  Soderberg 
correlation  only  accommodates  a  gas  deflection  of  up  to  140  . 
Since  the  efficiency  decreases  rapidly  at  angles  greater  than 
this,  no  harm  is  done  in  limiting  gas  deflection  to  140  . 

Reaction  is  another  parameter  which  is  limited.   The 
upper  bound  has  been  discussed  earlier  in  section  2.3.3.   A 
check  must  be  made,  however,  to  insure  that  at  no  point  does 
reaction  become  negative.   Physically,  a  negative  reaction 
would  mean  that  within  a  given  stage,  at  some  point,  a  pressure 
rise  was  occuring  to  be  compensated  for  by  a  larger  pressure 
drop  at  another  point.   The  overall  condition,  of  course, 
would  be  a  pressure  drop.   In  order  to  avoid  this  condition 
which  brings  about  higher  losses,  a  check  should  be  made  to 
insure  non-negative  reaction.   Negative  reaction  can  only 
occur  with  large  Vfi  and  since  it  is  assumed  that  r  V„  = 
constant,  the  check  should  be  made  at  the  smallest  radius, 
the  last  stage  hub. 

In  solving  for  a  mean  radius,  in  section  2.3.4,  no  con- 
tinuity check  was  made.   without  a  check  for  continuity  it 
would  be  possible,  using  the  model  described,  to  compute  a 
negative  machine  diameter.   In  order  to  avoid  this  and  to 
allow  for  bearings,  shafting,  and  blade  attachment,  a  check 
will  be  made  to  insure  that  the  hub  radius  calculated  is 
greater  than  50  percent  of  the  tip  radius. 
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The  one  dimensional  compressible  flow  equations  used 
throughout  also  pose  another  limitation.   The  fluid  flow 
throughout  the  machine  is  assumed  to  be  subsonic,  and 
the  design  has  proceeded  such  that  no  consideration  has  been 
given  to  the  shock  effects  caused  by  supersonic  flow.   To 
insure  subsonic  flow,  a  check  will  be  made  on  the  two  largest 
velocities  to  insure  a  Mach  number  less  than  one.   The 
velocities  are  the  last  stage  hub  velocity  at  rotor  inlet  and 
the  last  stage  tip  relative  velocity  at  rotor  outlet. 

In  addition  to  these  above  limitations,  there  are  physical 
limitations  imposed  by  the  computer  programming  itself.   These 
are  listed  in  the  appropriate  appendices. 
2.4   Radial  Turbine 
2.4.1   General 

The  term  radial  turbine  is  used  to  describe  a  wide  range 
of  turbine  configurations.   In  the  general  sense,  it  describes 
a  turbomachine  where  the  fluid  enters  the  rotor  with  no  appre- 
ciable axial  component  of  velocity.   Whether  the  fluid  flows 
radially  inward  or  radially  outward  and  whether  it  leaves  the 
rotor  in  a  radial  or  axial  direction  depends  upon  the  context 
of  the  discussion.   The  choice  made  here  to  limit  the  radial 
design  to  that  of  only  one  configuration  was  based  on  practi- 
cality.  The  90  degree  inward  flow  radial  turbine,  in  which 
the  fluid  enters  the  rotor  in  a  radial  but  leaves  in  an  axial 
direction,  is  the  most  common  type  of  radial  turbine  (7).   it 
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has  the  advantage  of  the  strongest  structural  strength,  also 
a  factor  which  can  be  important  at  high  rotational  speeds. 
Only  single  stage  turbines  will  be  considered  as  these  are 
the  most  common  also. 

The  methodology  of  designing  the  radial  turbine  can  be 
summarized  in  the  following  four  steps:  rotor  design,  nozzle 
design,  scroll  design  and  efficiency  correlation.   The 
station  designations  are  as  shown  in  Figure  6. 

2.4.2  Output  Required 

The  parameters  necessary  to  describe  the  turbine  of 
Figure  6  can  be  narrowed  down  to  the  following  list  for 
preliminary  design: 

1)  Nozzle  Inlet  and  Outlet  diameters 

2)  Nozzle  blade  height 

3)  Rotor  Inlet  diameter 

4)  Rotor  Outlet  diameters  at  hub  and  tip 

5)  Maximum  Scroll  Area 

6)  Unit  length 

7)  Specific  speed 

8)  Efficiency 

9)  Velocity  triangles  at 

A)  Nozzle  Inlet 

B)  Rotor  Inlet  and  Outlet. 

The  velocity  triangles  for  rotor  outlet  will  be  described 
at  hub,  mean,  and  tip  locations. 

2.4.3  Input 

The  parameters  necessary  to  describe  a  radial  machine  can 
be  categorized  in  the  same  manner  as  the  axial  parameters 
(section  2.3.3):  those  which  describe  the  mass  flow  rate,  those 
which  describe  the  work  accomplished  per  stuge,  and  those  which 
describe  the  performance  of  a  stage. 
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Radial  Turbine  Station  Designations 


Stator  — ^ 


^  Machine  Axis 


Figure  6 

As  it  has  been  previously  assumed,  that  the  radial  turbine 
design  is  for  a  single  staj^e  machine,  the  work  accomplished  per 
stage  is  the  same  as  the  total  work  to  be  accomplished  by  the 
cycle  parameters,  section  2.2.3. 

The  mass  flow  rate  for  this  machine  can  also  be  described 
in  a  manner  similar  to  that  for  the  axial  machine  and  which 
is  presented  in  section  2.2.2.   The  variables  necessary  for 
this  representation  are  the  mass  flow  rate  and  the  number  of 
machines  operating  in  parallel  along  the  fluid  flow  path. 

Performance  estimation  is  more  standardized  for  radial 
turbines  than  axial  turbines.   Through  the  use  of  a  non-dimen- 
sional specific  speed  parameter  it  is  possible  to  correlate 
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efficiency  with  rotational  speed,  volumetric  flow  rate,  and 
work  accomplished,  for  turbines  designed  with  the  same  basic 
assumptions  and  geometric  relationships.  (8)   There  are  several 
specific  speed  parameters  in  use  which  will  be  discussed  in  a 
later  section;  however,  for  the  immediate  purpose,  rotative 
speed  is  a  necessary  input.   One  last  input  remains  and  that 
relates  to  an  assumption  to  be  made  in  the  succeeding  section. 
Without  explaining  this  assumption  in  detail,  it  is  sufficient 
to  say  that  the  number  of  radial  blades  of  the  rotor  should 
be  included  as  an  input. 

The  input  for  the  radial  design  can  be  summarized  as 
follows : 

.  l)  Mass  flow  rate  of  the  cycle 

2)  Number  of  units  in  parallel 

3)  Shaft  revolutions  per  minute 

4)  Number  of  radial  rotor  blades 

5)  Cycle  conditions  of  section  2.2. 

2.4.4   Rotor  Calculations 

The  first  step  in  determining  the  rotor  parameters  is  to 
determine  the  conditions  at  rotor  inlet  (Station  3).   Borrow- 
ing an  analysis  from  the  field  of  centrifugal  compressor 
design,  it  can  be  shown  that  the  optimum  performance  is  ob- 
tained when  the  following  relationships  hold  at  rotor  inlet.  (8) 

?=1-  2.0 


e 


U3 


3=  1 


where  Z  is  the  number  of  radial  rotor  blades. 
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This  factor  (^)  is  often  called  the  slip  factor.   If  the 

velocity  at  rotor  exit  is  assumed  to  be  entirely  in  the  axial 

direction  for  reasons  similar  to  those  given  for  the  axial 

turbine  in  section  2.3.4,  then  the  slip  factor  may  be  used  in 

the  Euler  equation  to  solve  for  the  radius  at  rotor  inlet. 

Ah   =?U_2 
o    '  3 


eo 
where  U   =YTRPM  r  /30. 

Substituting  the  constants  into  this  gives: 

1 
r-  =  1510.8233  Uh  ) 2. 
i  RPM     Cl~¥ 

The  condition  necessary  to  avoid  negative  velocities  at  rotor 
inlet  (2)  is: 

m3   x   2  YP . 
u3   -  Z 

Combining  this  expression  with  that  of  the  slip  factor,  it  is 

possible  to  determine  the  optimum  inlet  flow  angle,  oc    . 

Tancc3  =  V03  =  V^/^ 

Vm3    VTS 

oC  =  ArcTan  (Z^g) . 
3  (2vr) 

It  is  now  possible  to  determine  the  complete  velocity  triangles 

at  rotor  inlet,  Figure  7. 

Using  the  above  relationships  and  the  definition  for  U, , 

V  ,  =  2KU, 
m3    £-   3 
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U      =  YTRPM   r., 
■*  30         J 


V        =   994.0819    l/ih    )2 
m3        — — t, o    1 

L  (    \     )5 

Y3    =    Vm3/cos    (c3) 


f03    =    V3    Sin    ^ 


/?3    =   ArcTan(V03-U3)/Vm3) 
W3   =   Vffl3/cos(/53). 


Radial    Turbine   Velocity    Triangles 


Rotor    Inlet 


Rotor  Outlet 


Figure  7 
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Using  the  compressible  flow  relationships,  the  flow  area 
at  inlet  is  obtained.   The  flow  area  and  annulus  area  are 
again  related  by  the  cosine  of  the  inlet  flow  angle  as  shown 
in  Figure  2.   The  annulus  area  in  this  case,  however,  consists 
of  the  circumference  times  the  blade  height.   From  this  the 
blade  height  can  be  calculated. 

At  rotor  outlet,  the  assumption  is  made,  as  stated  pre- 
viously, that  the  tangential  velocity  is  zero.   Watanabe  (9) 
states  that  there  is  a  definite  ratio  of  the  mean  outlet 
radius  to  inlet  radius  which  will  give  the  peak  efficiency. 
This  ratio  is  approximately  0.6.   It  has  also  been  somewhat 
traditional  to  design  for  a  relative  outlet  angle  at  the 
shroud  of  minus  60  .  (2)   This  represents  the  compromise  be- 
tween a  low  RPM  machine  and  a  low  loss  machine.   Applying  both 
of  these  criteria  in  a  modified  way,  an  outlet  velocity  tri- 
angle at  mean,  with  a  relative  flow  angle  of  minus  60   was 
chosen.   This  later  proved  to  correlate  well  with  previous 
work  as  shown  in  section  3.6  and  avoided  an  iterative  proce- 
dure for  calculating  rotor  outlet  conditions. 

Using  these  assumptions,  the  outlet  velocity  at  mean,  in 

terms  of  input  variables,  is 

1 

V  .  =  9.1347(Ah  A)2  (Radius  Ratio) 
m4  o  i 

The  value  of  0.6  specified  above  is  used  as  a  first  estimate 

for  radius  ratio.   If  the  resulting  outlet  conditions  do  not 

meet  the  checks  specified  in  section  2.4.9,  then  the  ratio  is 
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increased  and  outlet  conditions  recalculated. 

As  £.       is  minus  60  ,  the  other  outlet  conditions  at 

'  4m  ' 

mean  are 

IT   =VTRPM  r. 
4m    — q-   4m 

W  .   =  V.  /Cos  {/?.    )  . 
4m     4m'      Mm 

The  compressible  flow  relations  are  again  used  to  de- 
termine the  outlet  flow  area.   In  this  case,  since  the  outlet 
velocity  is  in  only  an  axial  direction,  the  flow  area  is 
equal  to  the  annulus  area.   The  outlet  height  may  be  calcula- 
ted by  the  same  method  used  in  section  2.3.5  to  calculate 
axial  blade  height: 

outlet  height  =  Area  ann/2y*r.  . 

Making  the  assumption  that  the  outlet  velocity  is  con- 
stant with  respect  to  radius,  and  knowing  the  outlet  hub  and 
tip  radii,  the  velocity  triangles  for  these  points  are  calcu- 
lated as  previously  shown,  with  the  exception  that  the  rela- 
tive outlet  angle  varies  with  the  radius,  or  more  simply 

/?  =  ATAN  (  -U7V.  ). 
4  4   4m 

2.4.5   Stator  Calculations 

The  stator  calculations  are  relatively  simple.   Watanabe 

(9)  states  that  the  clearance  between  stator  outlet  and  rotor 

inlet  should  be  twice  the  passage  height.   If  it  is  assumed 

that  the  blade  or  passage  height  is  the  same  for  both  rotor 

and  stator,  then  the  stator  outlet  radius  is  simply  the  radius 
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at  rotor  inlet  plus  two  times  the  blade  height.   As  in  most 

cases  this  gap  will  not  be  large,  and  as  it  can  be  assumed 

that  the  moment  of  momentum,  remains  constant  across  the  gap, 

r   =  V 
12    _03 

r3    V62 

where  r2/r^ « 1 . 

This  means  that  the  velocity  triangle  at  stator  outlet  is  sub- 
stantially the  same  as  at  rotor  inlet  and  will  not  be  recalcu- 
lated. 

The  radial  spacing  between  blades  or  pitch,  assuming  that 
there  are  as  many  stator  blades  as  rotor  blades,  is  the  circum- 
ference divided  by  the  number  of  blades.   The  throat  opening 
can  also  be  determined  in  a  manner  similiar  to  that  of 
section  2.3.5.   Using  a  pitch  to  chord  ratio  of  0.6,  the  blade 
chord  may  be  calculated.   Using  the  assumption  that  Cos^  - 
opening/spacing,  the  approximate  radius  at  the  stator  inlet 
can  be  calculated  as  shown  in  Figure  8. 

Angle  1  =  Yr  -  )r 
2     Z 

Angle  2  =  JT   -  ©C 
2 

Angle  3  =  \r  +  cC 

Z     b 

Side  1   =  r2  +  C  Cos  (>f  +  oC^ ) 

CI 

Side  2   =  C  Sin  {>T  +   cxl  ) 

Z     J 

radius  =    (Side  l2  +  Side  2~)~2 


Radial  Turbine  Stator  Configuration 
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Side  1 


Axis  of  Machine  into  page 


Figure  8 


The  inlet  angle  can  now  be  determined  using  the  law  of  Sines: 

Sin(c^)  =  r2  Sin(>r(  1-1  )-o| )  . 

'  Li 

ri 

For  the  passage  from  stator  inlet  to  outlet,  it  can  be  assumed 

that  the  fluid  is  incompressible.   This  is  a  valid  assumption 

as  the  Mach  number  is  low  at  rotor  inlet.   Since  the  mass  flow 

rate  must  remain  constant,  the  velocity  change  is  proportional 

to  the  flow  urea  change.   This  can  be  expressed  as 

V   =  V   r 
1     3   3  ' 
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Having  determined  the  inlet  angle, «C,  the  rest  of  the  velocity 

triangle  is  also  easily  obtainable. 

2.4.6   Scroll  Calculations 

The  scroll  shape  assumed  is  shown  in  Figure  9.   It  is 

assumed  that  the  area  is  distributed  radially  in  a  linear 

relationship,  thus: 

Area„  =  A     (1-6  ) 
8     max     ~T\r 

where  0  =  radial  angle  in  radians. 
The  maximum  area  needed  can  be  approximated  if  it  is 
assumed  that  the  tangential  and  meridional  components  of 
velocity  do  not  vary  as  a  function  of  the  radial  angle.   The 
maximum  area  is  then  equal  to  the  flow  area  at  stator  inlet 
which  in  turn  is  related  by  the  cosine  of  the  inlet  angle. 

Radial  Turbine  Scroll  Configuration 


:>Axis  of  Machine 


Figure  9 
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A     =2  r,YTt  cosk  ) 
max       1         1 

where  t  is  the  blade  height  at  rotor  inlet 

2.4.7   Radial  Efficiency 

As  with  the  axial  turbine,  the  theoretical  determination 
of  efficiency  for  a  radial  turbine  is  a  much  debated  subject. 
Wood  (4)  and  others  have  correlated  efficiency  with  the  non- 
dimensional  specific  speed  mentioned  previously,  section  2.4.3; 
however,  this  does  not  give  an  exact  enough  answer  for  the 
purposes  here.   A  more  analytical  method  is  needed. 

Heitt  &   Johnston  (10)  relate  the  performance  of  several 
specific  turbines  to  several  design  variables;  however,  the 
method  is  rather  a  narrow  one  intended  to  experimentally 
evaluate  performance  as  a  function  of  specific  variables, 
rather  than  provide  an  analytic  efficiency  determination. 
Rodgers  (11)  proposes  a  method  of  determining  efficiency 
which  accounts  for  losses  caused  by: 

1)  Disc  friction 

2)  Nozzle  shape 

3)  iiotor  friction 

4)  Rotor  blade  loading 

5)  Rotor  secondary  flow 

6)  liotor  incidence 

7)  Rotor  clearance 

8)  Rotor  cooling 

9)  Exhaust  diffuser. 

This  method,  while  verified  to  within  2  percent,  takes  into 
account  considerably  more  variables  than  are  necessary  for 
the  purposes  of  this  mathematical  model. 

The  method  chosen  was  that  used  by  Rohlik  ^8)  in  his 
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study  of  the  effects  of  turbine  geometry  on  efficiency.   This 
method  is  a  compilation  of  several  experimental  correlations 
relating  the  efficiency  to: 

1)  Stator  and  boundary  layer  losses 

2)  Clearance  losses 

3)  Rotor  windage  losses 

4)  Uotor  boundary  layer  losses 

5)  Exit  velocity  losses 

This  method  provided  the  needed  precision  and  was  based 
on  the  work  of  many  acknowledged  experts  in  the  field; 
Balje'  (12),  Heitt  k   Johnston  (11),  Wood  (7),  Stewart,  Vhitney 
&  Vang  (13)  and  Shepherd  (14)  to  name  a  few.   The  correction 
to  account  for  exit  velocity,  however,  is  omitted;  thus  the 
efficiency  calculated  is  a  total  to  total  efficiency  and  is 
compatible  with  the  efficiency  calculated  for  the  axial  tur- 
bine.  If  a  total  to  static  efficiency  is  desired,  it  can  be 
easily  calculated  from  the  data  given  for  the  outlet  velocity 
triangles  as: 

7  ts  =  jVt 

0    +  V   2 
<tt     4m 


2g  J(h  ,-h  0  ) 
6o     ol   o2s 

where  H ,  ,  =  efficiency  (total  to  total) 

').       =    efficiency  (total  to  static) 
x  s 

This  conversion  factor  can  also  be  applied  to  determine  total 

to  static  efficiency  for  the  axial  turbine  by  using  V    in 

lieu  of  V.  . 

4m 
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The  mechanics  of  the  method  consist  of  relating  each  loss 
to  an  associated  enthalpy  change  and  then  summing  these  in 
a  fractional  way  to  determine  efficiency.   The  first  loss  con- 
sidered is  that  of  rotor  windage.   Rotor  inlet  density  is  re- 
quired and  also  an  absolute  viscosity  Lt()  ,  at  rotor  inlet. 
The  density  is  computed  using  the  one-dimensional  compressible 
flow  relationships. 

/°  =  /o 

2 


/°   =  P  /RT 
o     o    o 


/>   =  P   (l-  Jf-1  M  ) 

RT       * 


2  i/(r-i)/  2    4  1/(^-1) 


For  Mach  numbers  less  than  one,  the  denominator  approa- 
ches one  and  the  density  is  then: 


,  i/(avi) 

2 


o 

RT 


Using  this  density,  the  absolute  viscosity,  Li()    can  be 
calculated  using  the  kinematic  viscosity  (S)  given  as  an 
input . 

A  Reynolds  number  for  the  rotor  inlet  is  defined  as: 

Substituting  in  the  constants  and  the  definition  of  rotative 
v  e  1  o  c  i  ty : 


55 


R  =  0.2094395  UPM  r  2/s  , 
e  3    3 

The  windage  loss  can  now  be  computed  as: 

2 


Lv  =  0.56/^iy 


R  *2  m.  106 
e     x 


or 


Lv  =  2.2^32^r3  < 


R  *2  m.  106 
e     i 

The  clearance  loss  is  expressed  as: 

L   =   h  VI 
c      o 

where  M.  is  defined  as 
1A   =  £l.004r,  +  .005r.,)  , 
t3        l4 
The  rotor  boundary  layer  loss  equations  involve  several 
assumptions  (8)  which  go  beyond  the  scope  of  this  model; 
however,  they  are  consistent  with  good  practice  and  those 
made  previously  in  this  model.   It  is  assumed  that  the  re- 
sulting equations  apply  to  this  model.   This  holds  true  for  the 
stator  boundary  layer  loss  equations  also. 

The  average  of  inlet  and  outlet  blade  spacing  is  defined 
as: 

S  ^(r3+r4h+r4t)- 

The  solidity  is  defined  as: 

fl*=  1.6r.   ((r.  /r.J-1) . 
4m     4  m   3 
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A  loss  coefficient  (e  )  is  calculated  by: 

e   =  (0.017  c-) (1  +  1. 9S   ) 

r    (0. 5-0. 003Z-0. 017a- )    rC+t  ) 

and  the  total  boundary  layer  loss  is  given  by: 

2 


L   =  e 
r     r 


\i 


4m 


1-e    2g   J 
r    fto 

In  a  like  manner  for  the  rotor,  a  stagger  angle  is  defined 

as : 

°st  =  <-  +*3)/2 

where  oc   is  given  by 

oc      =    ArcTan(  Sin^s. 
(- 


.  2t^  +  Co  s<Xo  \ 


)  ' 


The  loss  coefficient  (e  )  is  calculated  from 

s 


0.0076 


(1+Cos  c<  ) 
Si  , 


Cosc3-0.02  5  (   — oTT) 
and  a  total  loss  for  the  stator  boundary  layer  calculated 
from 


(V- 


) 


1-e   (2g   J) 
s    Bo 

The  efficiency  is  found  from 

)■).,     -Ah    -L  -L  A-^h  +L  +L  ). 
I XX  o   w   c'    o   s   r 

2.4.8   iNii  seel  laneous  Parameters 

Specific  speed  seems  to  be  defined  differently  by  each 

researcher.   Some  definitions  are  non-dimensional  and  others 
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are  not.  Wood  (7)  and  Rohlik  (8)  define  their  specific  speed 
similiar  to  that  used  in  centrifugal  pump  theory  as 

s       Vi:4/ 
H3'4 

where  Q,    is  the  volumetric  flow/sec 

H  is  the  ideal  head  drop  across 
the  turbine  in  feet. 

This,  however,  is  dimensional  and  for  that  reason  the 

non-dimensional  specific  speed  recommended  by  Carmichael  (2) 

is  used  instead. 

This  is  defined  as 

l 

N  =  RPM  Q,2 
s    60 


<«.  J*V3/4 

The  rotor  outlet  volumetric  flow  can  be  calculated  as 

^4  =  V/>4 

where  p.    is  obtained  using  the  compressible 
flow  relationship  of  the  previous  section. 

As  a  matter  of  interest  the  specific  speed  of  Wood  et  al 

may  be  obtained  from  that  used  in  the  model  by  using  the 

following  formula: 

N   (Vvood)  =  810.04455  N   (Uarmichael ) . 
s  s 

The  axial  unit  length  of  the  rotor  is  a  parameter  which 

must  be  calculated  in  detail  for  the  final  design;  however,  in 

order  to  simplify  the  calculations  in  the  mathematical  model 

a  simpler  method  is  used.   Based  on  the  data  given  in  several 
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of  the  references,  on  the  average,  the  unit  axial  length  is 
on  the  order  of  the  inlet  radius  less  the  outlet  hub  radius, 
thus  for  the  purposes  here, 

length  =  r3-r4h- 

The  machine  diameters  given  as  output  are  simply  calcu- 
lated from  the  previously  derived  radii. 
2.4.9   Radial  Limitations 

The  specific  speed  parameter  has  an  important  influence 
on  the  shape  of  the  turbine.   In  order  that  the  model  not 
generate  dimensions  which  were  unacceptable  from  a  geometric 
point  of  view  (i.e.  a  radical  unconventional  design  in  which 
other  factors  such  as  shroud  clearance  might  have  an  unaccount- 
ed for  effect),  the  input  parameters  are  first  checked  to  de- 
termine the  specific  speed  of  the  proposed  turbine.   If  this 
specific  speed  is  below  0.05  or  above  0.25  the  input  para- 
meters will  not  produce  a  design  which  is  conventionally 
feasible. 

The  calculation  of  stress  in  the  rotor  is  a  complicated 
subject  and  one  which  is  beyond  the  scope  of  the  preliminary 
design  phase.   Wood  (7),  however,  provides  a  means  of  deter- 
mining whether  a  design  will  be  stress  limited.   Specific 
speed  is  written  in  terms  of  outlet  area,  RPM,  rotor  tip  speed 
and  several  fluid  velocities.   Maximum  expected  values  for  a 
stress  limited  case  are  then  used  to  solve  for  a  specific 
speed.   This  specific  speed  was  approximate ly  twice  that 
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which  would  be  allowed  by  the  previous  criteria  and,  therefore, 
the  assumption  is  made  that  the  radial  turbines  allowed  by 
this  model  will  not  be  stress  limited. 

As  in  the  axial  turbine  a  continuity  check  must  be  made 
at  rotor  exit  to  assure  that  negative  diameters  are  not  calcu- 
lated.  The  rotor  outlet  radius  is  checked  to  insure  that  it 
is  at  least  10  percent  of  the  inlet  radius.   This  differs 
from  the  previous  50  percent  used  in  the  axial  turbine;  however, 
the  radial  rotor  doesn't  have  to  incorporate  a  through  shaft 
and  thus  can  be  of  a  much  smaller  inner  diameter.   If  the 
radius  in  question  is  too  small,  the  radius  ratio  used  in 
calculating  V   ,  section  2.4.4,  is  increased  by  2  percent. 
This  has  the  effect  of  increasing  the  density  of  the  fluid 
and  thus  decreasing  the  outlet  annulus  area  required. 

The  calculations  described  here  for  the  radial  turbine 
do  not  account  for  any  shock  effects  caused  by  operation  at 
a  Mach  number  greater  than  1.  luven   without  this  correlation, 
however,  Uodgers  (11)  states  that  Mach  number  effects  are 
not  important  until  M  =  1.2.   Carmichael  (2)  recommends  a 
limit  of  1.1.   The  lower  limit  was  used  to  be  consistent 
with  other  assumptions  made  previously.   The  two  velocities 
likely  to  be  the  largest,  the  outlet  tip  relative  velocity, 
and  the  rotor  inlet  velocity  are  checked  to  insure  that  they 
are  below  a  Mach  number  of  1.1.   If  not,  the  input  parameters 
result  in  an  unfeasible  turbine. 
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EVALUATION  OF  THE  MODEL 
3.1   General 

In  order  to  facilitate-  the  preliminary  design  of  turbines 
using  the  mathematical  model  developed  in  the  previous  chapter, 
it  was  necessary  first  to  program  the  model,  and  then  to  run 
the  program  over  a  wide  range  of  inputs  in  order  to  verify  the 
results  with  existing  correlations;  or  in  the  case  of  the  sim- 
pler geometric  relationships,  with  hand  calculations.   The 
areas  which  were  examined  for  validity  and  the  methods  used 
to  establish  same  can  be  summarized  as: 

Area  Method 

Simple  mathematical  hand  calculations 

relationships 

Aspect  ratio  selection  hand  calculations  of 

existing  turbine 

Soderberg  correlation  comparison  with 

function  original  function 

Axial  turbine  efficiency  comparison  with  other 

published  correlations 

Hadial  turbine  efficiency  comparison  with  other 

published  correlations 

Each  of  these  areas  will  be  described  in  more  detail  in  the 

following  sections. 

The  cycle  chosen  for  the  validation  was  the  Zener  sea 

cycle  as  proposed  by  Clarence  Zener  (15).   This  cycle  is 

similiar  to  others  (16)  which  seek  to  take  advantage  of  the 

temperature  differences  existing  in  the  world's  oceans,  and 

is  representative  of  the  type  of  cycle  for  which  the  raatlie- 
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matical  model  was  designed  to  be  used  to  evaluate.   The  de- 
tails of  this  cycle  and  of  the  individual  parameters  which 
were  examined  are  contained,  in  Appendix  F.       It  is  sufficient 
to  say  here  that  over  50  different  combinations  of  input  para- 
meters were  examined  for  each  of  two  different  working  fluids. 
This  range  of  input  variables  provided  a  representative 
picture  of  the  design  surface,  and  allowed  the  results  to  be 
broad  enough  to  correlate  with  existing  data  over  a  wide 
range . 

3.2  Simple  Mathematical  Relationships 

This  facet  of  the  model  is  the  predominent  one,  as 
explained  in  section  1.3,  and  also  the  easiest  to  verify. 
The  major  portion  of  the  model  is  based  on  traditional  theory 
which  has  been  well  proven  in  the  laboratory.   It  is  not  the 
purpose  here  to  attempt  to  revalidate  the  conclusions  of  re- 
searchers far  more  qualified  than  the  author.   The  emphasis 
has  been  placed  instead  on  choosing  accepted  relationships 
and  correlating  them  into  a  usable  model. 

Based  on  this,  the  validation  of  this  segment  consisted 
of  hand  calculating  the  various  mathematical  relationships 
and  verifying  the  results  of  the  computer  program.   To  this 
end,  the  model  has  been  debugged  of  all  forseeable  mechanical 
type  mathematical  errors. 

3.3  Aspect  hatio  Correlation 

The  decision  to  develop  an  internal  choice  of  aspect 
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ratio  based  on  blade  bending  stress  for  the  axial  turbine, 
rather  than  have  the  designer  predetermine  it,  was  based  on 
primary  and  secondary  input,  factors  as  explained  in  section 
2.3.3.   It  was  necessary  then  to  verify  that  the  internal 
choice  was  a  reasonable  one,  consistent  with  existing  designs. 

This  was  carried  out  by  applying  the  correlation  developed 
in  section  2.3.4  to  an  existing  design  for  an  intermediate 
stage  turbine,  developed  for  marine  propulsion  use  and  des- 
cribed by  T.  B.  Hutchinson  in  reference  17.   This  design  was 
for  an  8  stage  steam  turbine  operating  at  5036  RPM  and  pro- 
ducing 8940  horsepower  at  90  percent  efficiency  with  an  inlet 
temperature  and  pressure  of  1000  F  and  225  psig  respectively. 
Using  the  data  supplied  by  Hutchinson  in  reference  17,  and 
making  assumptions,  where  necessary,  of  input  variables  which 
would  optimize  efficiency,  the  correlation  of  section  2.3.4 
was  applied  to  the  design  and  an  aspect  ratio  was  computed. 
This  aspect  ratio  was  approximately  1.4.   The  aspect  ratio 
chosen  in  the  actual  design  was  approximately  1.5. 

While  this  simple  point  verification  does  not  guarantee 
validity  over  the  complete  range  of  input  variables,  it  does 
show  that  the  relationship  is  reasonable.   When  this  relation- 
ship is  combined  with  an  upper  limit  for  aspect  ratio,  to  take 
into  account  the  fact  that  other  factors  than  stress  may  limit 
the  design,  the  resulting  range  of  variation  of  aspect  ratio 
(0.25  to  10.0)  is  consistent  with  existing  designs  (5)  and  is 
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assumed  to  provide  a  reasonable  result. 

3.4  Soderberg  Correlation 

In  determining  efficiency  for  the  axial  turbine,  section 
2.3.5,  it  was  necessary  to  formulate  a  function  which  approx- 
imated the  curve  of  emperical  data  given  by  Soderberg  (4)  for 
the  correlation  of  gas  deflection  angle  and  the  primary  loss 
factor.   This  was  accomplished  in  a  computer  subprogram 
through  the  use  of  LaUrange  polynomials  to  curve  fit  the 
given  data  to  a  third  order  approximation  of  the  function 
(18).   Without  explaining  the  details  of  the  theory  of  curve 
fitting,  it  is  sufficient  to  state  that  within  a  given  sub- 
interval  the  correlation  was  represented  by  the  function 
W  =  X  Y,  ■+  X2  Y2  +  X   Y 

where  X,  ,  X~  and  X-,  are  third  order  functions  dependent  upon 
the  value  of  the  abscissa  given?  and  Y,  ,  Y„,  and  Y_.  are  values 
of  the  ordinate  which  bracket  the  given  abscissa  value. 

The  simulation  of  this  function  is  shown  in  Figure  10. 

3.5  Axial  Turbine  Efficiency 

Although  the  correlation  adopted  for  determination  of 
axial  total  to  total  efficiency  has  been  previously  validated, 
it  was  necessary  to  compare  the  results  obtained  with  the 
model  to  other  correlations  in  order  to  assure  that  the 
correlations  provided  reasonable  results  for  the  assumptions 
made  in  the  model.   This  was  accomplished  by  comparing  the 
computed  efficiency  of  several  example  cases  for  the  Zener 
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cycle  with  the  correlations  presented  by  Craig  and  Cox  (3). 
This  comparison  is  shown  in  Figure  11. 

This  comparison  does  not  account  for  the  differences 
between  the  two  correlations  ixi  their  accounting  for  secon- 
dary loss  parameters;  however,  the  results  are  close  enough 
to  show  that  the  application  of  Soderberg's  method  is  valid 
and  consistent  with  accepted  practices. 
3.6   Radial  Turbine  Efficiency 

It  was  also  necessary  to  compare  the  computed  radial 
turbine  efficiency  to  other  correlations  in  order  to  verify 
its  application  in  the  model.   This  can  be  more  easily  accom- 
plished in  the  radial  turbine  than  the  axial  through  the  use 
of  the  non-dimensional  specific  speed  parameter  mentioned 
in  section  3.4.8.   Several  examples,  covering  a  wide  range 
of  specific  speeds,  were  compared  with  the  correlations  of 
Wood  (2)  and  Rohlik  (8).   The  results  of  this  correlation 
are  shown  in  Figure  12.   The  close  agreement  over  the  wide 
range  of  specific  speeds  shows  that  the  determination  of 
radial  efficiency  is  also  valid  and  consistent  with  accepted 
practice . 


Axial  Turbine  Efficiency  Comparison 
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Craig  &  Cox  Correlation 
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Radial  Turbine  Efficiency  Comparison 
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CONCLUSIONS 

The  formulation  of  conclusions  about  a  project  of  this 
nature  is  a  somewhat  limited  undertaking.   The  object  of  the 
thesis  is  to  compile  existing  data,  relationships,  and  methods 
in  the  field  of  turbine  design  and  to  formulate  and  compute- 
rize a  mathematical  model,  based  on  this  compilation.   This 
has  been  accomplished  as  shown  in  chapter  2. 

The  validity  of  the  mathematical  model  has  been  estab- 
lished in  chapter  3  to  the  extent  that  it  is  valid,  consistent 
with  the  assumptions  stated  in  its  formulation,  and  within 
the  requirements  posed  by  the  preliminary  design  stage. 

The  actual  use  of  this  programmed  model,  as  explained 
in  Appendix  F,  has  shown  that  through  the  utilization  of 
this  program,  or  programs  of  this  type,  the  designer  may 
explore  a  wide  range  of  input  variables  in  a  comparably 
short  amount  of  time.   This  enables  him  to  more  fully  examine 
the  design  spectrum  and  to  develop  a  qualitatively  better 
estimation  of  thermal  cycle  feasibility. 
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HECOMMENDATJONS 

In  developing  the  mathematical  model  described  in  chapter 
2,  it  became  apparent  that  there  were  four  areas  where  the 
model  could  be  expanded  which  would  increase  its  usefulness 
to  the  designer.   The  first  of  these  is  in  the  area  of  the 
flow  relationships  utilized.   The  one  dimensional  flow  rela- 
tionships do  not  completely  describe  the  flow.   Expansion  of 
these  relationships  to  include  shock  effects  would  broaden 
the  applicability  of  the  model. 

The  second  and  third  areas  which  could  be  expanded  are 
integrally  related  to  one  another.   They  entail  enlarging 
the  areas  of  the  machine  which  are  described  by  the  model, 
perhaps  including  material  limits  and  selection,  disc  design, 
bearing  and  shafting  design.   Concurrent  with  this  would  be 
the  increasing  number  of  input  variables  necessary  to  des- 
cribe the  machine.   This  would  be  limited  by  the  amount  of 
detail  required  to  satisfy  a  particular  designer's  definition 
of  preliminary  design. 

The  last  area  which  was  not  included  in  this  model  is 
that  of  cost  estimation.   Information  in  this  area  is  gene- 
rally proprietory  and  not  easily  obtainable.   host  of  the 
information  available  is  limited  in  scope  and  not  applicable 
to  the  general  case  as  described  herein.   This  area  in  itself 
could  be  the  subject  of  a  dissertation,  but  certainly  more 
work  is  needed  before  it  is  acceptable  for  inclusion  in  a 
model  similar  to  that  developed  in  chapter  2. 
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APPENDIX  A 

A.l   Variable  List 

The  following  is  a  list  of  variables  as  they  appear  in 
the  main  program  and  the  subprograms.   They  are  listed  in  the 
order  in  which  they  appear  within  the  body  of  the  programs. 
Subscripted  variables  refer  to  machine  inlet  and  outlet, 
subscripts  1  and  2  respectively.   The  units  of  these  variables 
are  those  listed  in  the  Notation  section,  pages  7  to  9,  with 

the  exception  of  pressure,  kinematic  viscosity  and  all  angles 

2     2 
which  are  internally  converted  to  lbf/ft  ,  ft  /sec  and 

radians  respectively. 

MAIN  PROGRAM 

Cycle  and  Machine  Inputs 

HOI  inlet  stagnation  enthalpy 

TO  stagnation  temperature 

PO  stagnation  pressure 

NU  kinematic  viscosity 

1102  outlet  stagnation  enthalpy 

H02S  outlet  isentropic  stagnation  enthalpy 

T02S  outlet  isentropic  stagnation  temperature 

P02S  outlet  isentropic  stagnation  pressure 

OP  specific  heat  at  constant  pressure 

GAMA  ratio  of  specific  heats 

0V  specific  heat  at  constant  volume 
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H  universal  gas  constant 

DHO  stagnation  enthalpy  drop  across  machine 

HJ  universal  gas  constant  times  778 

MM  number  of  combinations  of  input  parameters  entered 

AXN  number  of  turbine  stages-axial 

REAM  reaction  at  mean 

RPMX  shaft  revolutions  per  minute-axial 

ALFAMN  rotor  inlet  true  flow  angle  at  mean 

UN1TX  number  of  axial  turbines  in  parallel  on  flow  path 

RPMR  shaft  revolutions  per  minute-radial 

UN1TR  number  of  radial  turbines  in  parallel  on  flow  path 

ZR  number  of  radial  blades 

WW  total  mass  flow  rate  of  working  fluid 

Axial  Turbine 


ft-ADM        radius  at  mean 

VTMN        tangential  velocity  component  at  rotor  inlet, 
mean  radius 

ViUN  true  velocity  at  rotor  inlet,  mean  radius 

VX  axial  velocity  component 

WMR  relative  velocity  at  rotor  outlet,  mean  radius 

TALLIN  tangent  of  ALJ7AMN 

BETAMN  relative  flow  angle  at  rotor  inlet,  mean  radius 

BETAhil  relative  flow  angle  at  rotor  outlet,  mean  radius 

ZETANl      primary  boderberg  efficiency  parameter  at  rotor 
inlet 


74 


ZETARi      primary  Soderberg  efficiency  parameter  at  rotor 
outlet 

CON         product  or  radius  times  tangential  velocity 
component 

WMN  relative  velocity  at  rotor  inlet,  mean  radius 

UMN  tangential  rotor  velocity  at  mean  radius 

FLFN  function  f  (section  2.2.3) 

AFLO  flow  area 

AANN  annulus  area 

T  blade  height-axial 

KADT  tip  radius 

RADII  hub    radius 

VTTN        tangential  velocity  component  at  rotor  inlet, 
tip  radius 

VTHN       tangential  velocity  component  at  rotor  inlet, 
hub  radius 

ALFATN  rotor  inlet  true  flow  angle  at  tip  radius 

ALFAHN  rotor  inlet  true  flov  angle  at  hub  radius 

VTN  true  velocity  at  rotor  inlet,  tip  radius 

V11N  true  velocity  at  rotor  inlet,  hub  radius 

UTN  tangential  rotor  velocity  at  tip  radius 

UHN  tangential  rotor  velocity  at  hub  radius 

WTTN        relative  tangential  velocity  component  at  rotor 
inlet,  tip  radius 

WTI1N        relative  tangential  velocity  component  at  rotor 
inlet,  hub  radius 

BKTATN      relative  flow  angle  at  rotor  inlet,  tip  radius 

BETAHN      relative  flow  angle  at  rotor  inlet,  hub  radius 
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wTN  relative  velocity  at  rotor  inlet,  tip  radius 

WHN  relative  velocity  at  rotor  inlet,  hub  radius 

BETATR  relative  flow  angle  at  rotor  outlet,  tip  radius 

BETAHR  relative  flov  angle  at  rotor  outlet,  hub  radius 

WTH  relative  velocity  at  rotor  outlet,  tip  radius 

Willi  relative  velocity  at  rotor  outlet,  hub  radius 

CHD  blade  chord 

ASltAT  aspect  ratio 

SPN  stator  spacing 

TBTAMN  tangent  of  BETAMN 

TBTAMR  tangent  of  BETA^iR 

SPR  rotor  spacing 

TR  throat  opening  for  the  rotor 

TN  throat  opening  for  the  stator 

DR  hydraulic  diameter  for  rotor 

DN  hydraulic  diameter  for  stator 

RER  Reynolds  number  for  rotor 

REN  Reynolds  number  for  stator 

ZETAR2  Soderberg  efficiency  parameter  corrected  for 
aspect  ratio  for  rotor 

ZETAN2  Soderberg  efficiency  parameter  corrected  for 
aspect  ratio  for  stator 

ZETAR3  Soderberg  efficiency  parameter  corrected  for  All 
and  Re,  for  rotor 

ZETAN3  Soderberg  efficiency  parameter  corrected  for  AR 
and  Re,  for  stator 


EFJ?1 


stage    efficiency 
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ZX 
NZX 
EFEX 
XMN 

XMR 

HEAH2 
ESTMX 
DIAM1 
DIAM2 
UNLEN 
u 


number  blades-decimal  number 

number  blades-whole  number 

machine  efficiency 

true  Mach  number  at  rotor  inlet,  hub  radius,  last 
stage 

relative  Mach  number  at  rotor  outlet,  tip  radius, 
last  stage 

reaction  at  last  stage,  hub  radius 

centrifugal  stress  at  blade  root,  last  stage 

inlet  diameter 

outlet  diameter 

machine  length 

radians  to  degrees  conversion 


itadial  Turbine 


Z 

KAU3 

ALFA  3 

VM3 

V3 

XM3 

FUN  3 

AFLU3 

IK 

KAU2 

itUAT 


slip  factor 

radius  at  rotor  inlet 

true  inlet  flow  angle,  rotor  inlet 

meridional  velocity  component  at  rotor  inlet 

velocity  at  rotor  inlet 

Mach  number  at  rotor  inlet 

function  f  at  rotor  inlet  (section  2.2.3) 

flow  area  at  rotor  inlet 

stator  blade  height 

radius  at  stator  exit 

radius  ratio  (section  2.4.4) 
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V4M  velocity  at  rotor  outlet,  mean  radius 

V4H  velocity  at  rotor  outlet,  hub  radius 

H.AD4M  mean  radius  at  rotor  outlet 

XM4M  Mach  number  at  rotor  outlet 

PUN4  function  f  at  rotor  outlet  (section  2.2.3) 

AANN4  annulus  area  at  rotor  outlet 

T4  blade  height  at  rotor  outlet 

RAD4H  hub  radius  at  rotor  outlet 

RAD4T  tip  radius  at  rotor  outlet 

BB  exponential  constant 

BEN4  fluid  static  density  at  rotor  outlet 

y  volumetric  flow  rate  at  rotor  outlet 

SPSP  non-dimensional  specific  speed 

SPRR  stator  outlet  spacing 

OP  stator  throat  opening 

UR  stator  chord 

RAD1  stator  inlet  radius 

SALPA1  sine  of  ALPA1 

CALPA1  cosine  of  ALPA1 

ASCR  maximum  scroll  area 

VI  velocity  at  stator  inlet 

V>il  meridional  velocity  component  at  stator  inlet 

VT1  tangential 

ALPA1  true  flow  angle  at  stator  inlet 

VT3  tangential  velocity  component  at  rotor  inlet 
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U3  tangential  rotor  velocity  at  inlet 

BETA3  relative  flow  angle  at  rotor  inlet 

W3  relative  velocity  at  rotor  inlet 

U4M  tangential  rotor  velocity  at  outlet,  mean  radius 

U4H  tangential  rotor  velocity  at  outlet,  hub  radius 

U4T  tangential  rotor  velocity  at  outlet,  tip  radius 

BETA4M  relative  flow  angle  at  rotor  outlet,  mean  radius 

BETA4H  relative  flow  angle  at  rotor  outlet,  hub  radius 

BETA4T  relative  flow  angle  at  rotor  outlet,  tip  radius 

W4M  relative  velocity  at  rotor  outlet,  mean  radius 

V4H  relative  velocity  at  rotor  outlet,  hub  radius 

W4T  relative  velocity  at  rotor  outlet,  tip  radius 

XiM4  relative  Mach  number  at  rotor  outlet,  tip  radius 

UE3  Reynolds  number  at  rotor  inlet 

DEN3  fluid  static  density  at  rotor  inlet 

XLM  windage  loss 

XLC  clearance  loss 

SPltM  mean  rotor  spacing 

SKiMAR  solidity 

Eli  rotor  boundary  layer  loss  coefficient 

XLK  rotor  boundary  layer  loss 

ALFAO  stator  angle  (section  2.4.7) 

ALFAST  stator  stagger  angle 

ES  stator  boundary  layer  loss  coefficient 

XLS  stator  boundary  layer  loss 
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EFFR  machine  efficiency 

UNLENR  machine  length 

DIA3  diameter  at  rotor  inlet 

DIA4T  tip  diameter  at  rotor  outlet 

UIA4M  hub  diameter  at  rotor  outlet 

JD1A1  diameter  at  stator  inlet 

D1A4M  mean  diameter  at  rotor  outlet 

DIA2  diameter  at  stator  outlet 

SUBPROGRAM  SODCOR 

XI  gas  deflection  angle,  radians 

X  gas  deflection  angle,  degrees 

XX  abscissa  for  Soderberg  correlation 

YY  ordinate 

XL1  first  LaGrange  Polynomial 

XL2  second  LaCirange  Polynomial 

XL3  third  LaGrange  Polynomial 

Y  Soderberg  gas  deflection  loss  parameter 

SUBPROGRAM  FLOREL 

GAMA   .  ratio  of  specific  heats 

R  universal  gas  constant 

UP  specific  heat  at  constant  pressure 

TO  stagnation  temperature 

V  velocity 
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Xh  Mach  number 

FLFN  function  f  (section  2.2.3) 

T  static  temperature 

EXPR  internal  constant 

CON  internal  constant 

FUN  internal  constant 
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APPENDIX  B 
B.l   ivlain  Program 

A  complete  listing  of  .the  main  program  is  contained  on 
the  following  pages.   Appendix  E    contains  detailed  instruc- 
tions for  input-output  and  running  the  program.   The  mathe- 
matical model  for  the  program  is  developed  in  chapter  2. 
The  logic  sequence  used  in  programming  the  model  is  outlined 
below : 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 

11 
12 

13 

14 
15 
16 

17 
18 


19 
20 
21 


Definition  of  real  variables  and  common  variables 
Dimensioning  of  subscripted  variables 
Input  format  specification 
Output  format  specification 
Readin  input  data 
Readout  input  data  for  checking 
Convert  input  data  to  standard  units 
Calculate  fluid  properties 
Readout  fluid  properties  for  checking 
Commence  do-loop  for  calculations  on  each  set  of 
turbine  input  parameters 

Calculate  axial  turbine  mean  parameters 
Calculate  axial  first  and  last  stage  parameters  on 
internal  do— loop 

If  design  is  rejected  on  specified  checks,  axial  cal- 
culations are  dropped  and  radial  are  commenced. 
Calculate  axial  machine  parameters 
Convert  axial  variables  to  output  units 
Readout  output  data  for  axial  machine 
Calculate  radial  parameters 

If  design  is  rejected  on  specified  checks,  radial 
calculations  are  dropped  and  next  set  of  axial 
calculations  commenced 

Convert  radial  variables  to  output  units 
iteadout  output  data  for  radial  machine 
Commence  calculations  for  next  set  of  axial  input 
parameters 
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APPENDIX  C 
C   Subprogram  SODCOR 

This  program  computes  an  efficiency  parameter  from 
given  data  points  using  LaGrange  Polynomials  in  a  curve- 
fitting  logic.   A  complete  listing  of  the  subprogram  is 
contained  on  the  following  page.   The  mathematical  model 
is  explained  in  section  3.4.   The  logic  sequence  used  in 
programming  the  model  is  outlined  below: 

1)  Transfer  of  input  variables 

2)  Definition  of  curve  data  points 

3)  Conversion  of  input  data  to  standard  units 

4)  Classification  of  input  data  as  to  sub— division 
of  curve  for  curve-fitting  purposes 

5)  Definition  of  sub-division  curve-fitting  data 

6)  Calculation  of  LaGrange  Polynomials 

7)  Calculation  of  output  variable 

8)  Transfer  of  output  variable 
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APPENDIX  D 
D   Subprogram  FLOHEL 

This  subprogram  computes  a  flow  function  and  Mach  number 
using  the  one-dimensional  flow  functions.   The  mathematical 
model  is  explained  in  section  2.2.3.   A  complete  listing  of 
the  program  is  contained  on  the  following  page.   The  logic 
sequence  used  in  programming  the  model  is  outlined  below. 

1)  Transfer  of  input  variables 

2)  Transfer  of  common  variables 

3)  Calculation  of  hach  number 

4)  Calculation  of  flow  function 

5)  Transfer  of  output  variables 
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APPENDIX  E 
E.l   Specific  Instructions  for  Use 

This  section  is  intended  to  provide  detailed  instruc- 
tions for  the  use  of  the  computer  programs  described  in  the 
previous  appendices  including  input-output  interpretation. 
The  specific  program  control  cards,  described  are  required 
on  the  Interdata  70  computer.   Check  the  user's  manual  be- 
fore running  on  other  machines  as  these  control  cards  may 
differ. 
E.2   Loading  Sequence 

The  sequence  for  loading  the  main  and  subprograms  with 
their  required  control  cards  is  as  follows: 

1)  Job  card 

2)  Language  card 

3)  Listing  designation  card 

4)  Main  program 

5)  Subprogram  SUDCOR 

6)  Subprogram  FLOREL 

7)  Execution  card 

8)  Common  length  designation  card 

9)  Data  cards 

10)   Termination  card 

E.3   Job  Control  Cards 

-    The  JCB  card,  language  card,  listing  designation  card 
and  JOB  termination  cards  are  all  of  the  standard  format, 
common  to  any  job  run  on  the  Interdata  70  system.   The  ex- 
ecution card  must  carry  a  numeral  one  in  column  seventeen  in 
addition  to  the  standard  format.   The  common  length  designa- 
tion card  specifies  the  space  in  blank  common  to  be  used  in 
the  job.   The  format  for  this  card  is  *BC  000C  commencing 
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in  column  one. 
E.4   Data  Cards 
£.4.1   General 

The  Data  cards  necessary  for  use  of  the  program  can  be 
grouped  into  three  categories. 

1)  Cycle  Data 

2)  Number  of  Design  Sets 

3)  Turbine  Design  input  sets 

They  must  be  loaded  in  the  above  order.   Card  groups  1  and  2 
are  only  required  once  per  job,  while  card  group  3  is  re- 
quired for  each  individual  design  set  examined. 
E.4. 2   Cycle  Data 

Three  cards  are  required  to  specify  the  data  for  a 
given  thermal  cycle.   The  context  of  each  is  given  below. 
No  information  can  bo  omitted,  and  they  must  be  loaded  in 
the  order  described. 
SPACES  INFORMATION  UNITS  FORMAT 


CARD  ONE 


1-: 

LO 

20- 

-30 

40- 

-50 

60- 

-70 

CARD  TWO 

TU/lbm 

F10.5 

°R 

F10.5 

psia 

F10.5 

ft2/hr 

F10.5 

inlet  total  enthalpy 
inlet  total  temperature 
inlet  total  pressure 
inlet  kinematic  viscosity 


1-10         outlet  total  enthalpy       BTU/lbm  F10.5 

20-30        outlet  total  temperature       R  F10.5 

40-50        outlet  total  pressure        psia  F10.5 

2 

60-70        outlet  kinematic  viscosity   ft  /hr  F10.5 


SPACES 
CARD  THREE 
1-10 

20-30 

40-50 


INFORMATION  UNITS 


outlet  isentropic  total    BTU/lbm 
enthalpy- 
outlet  isentropic  total       R 
temperature 

outlet  isentropic  total      psia 
pressure 
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FORMAT 

F10.5 
F10.5 
F10.5 


E.4.3   Number  of  Design  Sets 

One  card  is  necessary  to  specify  the  number  of  sets  of 
turbine  parameters  entered.   A  set  consists  of  one  combina- 
tion of  axial  turbine  input  parameters,  and  one  combination 
of  radial  turbine  input  parameters,  with  both  combinations 
sharing  the  same  total  mass  flow  rate  (section  2.2.2).   The 
program  will  accommodate  up  to  99  separate  sets  of  input 
data.   The  number  of  sets  must  be  specified  in  12  format  on 
spaces  one  and  two  of  the  card  with  the  other  spaces  left 
blank.   This  card  may  not  be  omitted. 
E.4.4   Turbine  Design  Input  Sets 

A  turbine  design  input  set  consists  of  a  combination  of 
axial  input  parameters,  section  2.3.3,  and  a  combination  of 
radial  input  parameters,  section  2.4.3,  which  share  the  same 
total  mass  flow  rate,  section  2.2.2.   Three  cards  are 
necessary  for  each  set  and  no  ini'ormati  on  or  card  may  be 
omitted.   They  must  be  loaded  in  the  order  in  which  they  are 
described.   In  the  event  that  a  run  is  not  to  be  designated 
by  name,  a  blank  card  must  be  included  in  place  of  card  one. 


F5.0 

F5.2 

Rev/rain 

•  F10.0 

Degrees 

F5.2 

*r 


100 

SPACES  INFORMATION  UNITS  FORMAT 

» 

CARD  ONE 

6-70  Any  desired  run  designa—  H 

tion  information 

CARD  TWO-Axial  Turbine 

11-15         Number  of  stages 

21-25         Inlet  reaction  at  midblade 

31-40         Shaft  RPM 

51-55         Rotor  inlet  flow  angle 
at  midblade 

61-65         Number  of  machines  op-  F5.0 

erating  in  parallel  on 
total  flow  path  (may  be  one) 

CARD  THREE-Radial  Inflow  Turbine 

11-20         Shaft  RPM  Rev/mi n        F10.0 

31-35         Number  of  machines  op-  F5.0 

erating  in  parallel  on 
total  flow  path  (may  be  one) 

41-45         Number  of  rotor  blades  F5.0 

51-60         Total  mass  flow  rate  for     lbra/sec        F10.2 
all  machines-i . e .  sum  of 
all  parallel  flow  paths. 
See  section  2.2.2.   It 
must  be  the  same  for  Axial 
and  Radial  Design  Set. 

E.5   Output 

E.  5 . 1   General 

The  cycle  data  will  be  printed  out  once  for  each  job, 

in  the  same  format  it  was  entered  in  with,  as  a  check  for 

accuracy.   The  universal  gas  constant,  in  units  of  ft  lbf/ 

lbm  R,  the  specific  heat  at  constant  pressure,  and  the  ratio 

of  specific  heats  will  also  be  printed  out  as  a  check.   For 
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each  turbine  design  data  set,  the  run  designation  information, 

# 

the  assumed  input  variables,  and  the  preliminary  design  data 
will  be  printed  out.   If  the  assumed  input  variables  result 
in  a  design  which  doesn't  pass  the  axial  or  radial  feasibi- 
lity checks,  sections  2.3.8  and  2.4.9,  the  design  data  will 
be  omitted  and  an  unf easibili ty  message  printed,  section 
E.5.3. 
E.5.2   Preliminary  Design  Data 

The  format  of  the  preliminary  design  data  is  self- 
explanatory.   The  output  is  in  units  of  feet  for  all  dimen- 
sions, degrees  for  all  angles,  and  feet  per  second  for  all 
velocities.   Stress  is  given  in  psi  and  tpi .   The  abbrevia- 
tions used  in  describing  the  velocity  triangles  are  as 
follows : 
Axial  Turbine 

UHN,  UivlN,  UTN,  -  rotor  tangential  velocity  at  hub,  mean  and 

tip  radii 

VTHN,VTMN,VTTN,-  fluid  tangential  velocity  component  at 

rotor  inlet  at  hub,  mean  and  tip  radii 

VX  -  fluid  axial  velocity  component 

VHN,  VMN,  VTN,  -  fluid  velocity  at  rotor  inlet  at  hub,  mean 

and  tip  radii 

WHN,  WMN,  VTN,  -  fluid  relative  velocity  at  rotor  inlet  at 

hub,  mean  and  tip  radii 

Willi,    MMR,    WTli,    -    fluid    relative    velocity    at    rotor    outlet,    hub, 

mean    and    tip    radii 

ALPAHN,ALFAMN,ALJJ'ATN    -    rotor    inlet    flow    angle    at    hub,    mean    and 

tip    radi  i 
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BETAHN,BETAMN,BETATN,  -  relative  rotor  inlet  flow  angle  at 

hub,  mean  and  tip  radii 

BETAHlt,BETAMK,BETATR,  -  relative  rotor  outlet  flow  angle  at 

hub,  mean  and  tip  radii 

Radial  Turbine 

VI,  V3         --  fluid  velocity  at  station  1  or  3 

VT1 ,  VT3,      -  fluid  tangential  velocity  component  at 

station  1  or  3 

VM1 ,  VM3 ,      -  fluid  meridional  velocity  component  at 

station  1  or  3 

U3  -  rotor  tangential  velocity  at  inlet 

W3  -  fluid  relative  velocity  at  rotor  inlet 

ALFA1 , ALEA3 ,  -  inlet  flow  angle  at  station  1  or  3 

BETA3  -  relative  rotor  inlet  flow  angle 

U4H,U4M,U4T,   -  rotor  outlet  tangential  velocity  at  hub, 

mean  and  tip  radii 

V4H,V4M,V4T,   -  rotor  outlet  fluid  velocity  at  hub,  mean 

and  tip  radii 

W4H,W4M,W4T,   -  rotor  outlet  fluid  relative  velocity  at 

hub,  mean  and  tip  radii 

BETA4H,BETA4M,BETA4T,  -  rotor  outlet  fluid  relative  flow 

angle  at  hub,  mean  and  tip  radii 

E.5.3   Unf easibili ty  Conditions 

If  a  given  set  of  turbine  design  inputs  does  not  meet 

the  feasibility  checks  described  in  sections  2.3.H  and  2.4.9, 

a  design  unfeasibile  message  will  be  printed.   As  a  reference, 

the  checks  associated  with  each  message  are  summarized  below. 

MESSAGE  INTERPRETATION 

Axi  al  Turb i  ne 

Gas  Deflection  too  big         Either  **   or  ^/?  -;-/£  )  is  greater 

than  140° 
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MESSAGE 

Hub  Radius  too  small 

Hub  iteal  Mach  Number 
too  big 

Tip  Real  Mach -Number 
too  big 

Hub  reaction  too  small 

Radial  Turbine 

Specific  Speed  to  Radical 


Inlet  Mach  number  too 
high 

Outlet  Relative  Mach 
number  too  high 


INTERPRETATION 


r,  <  0.5  r, 
h        t 


Last  stage  hub  true  Mach 
number  greater  than  0.99 

Last  stage  tip  relative  Mach 
number  greater  than  0.99 

Last  stage  hub  reaction  less 
than  0.01 


Machine  specific  speed  less 
than  0.05  or  greater  than 
0.25. 

Mach  number  at  rotor  inlet 
greater  than  1.1. 

Relative  Mach  number  at  rotor 
tip  greater  than  1.1 


E.5.4   Compilation  Errors 

It  is  possible,  given  certain  sets  of  axial  turbine  in- 
put data, to  cause  a  computer  run  time  error  (negative  square 
root  and  overflow).   This  condition  sometimes  occurs  during 
the  calculation  of  aspect  ratio.   It  has  been  forseen  in  pro- 
gramming the  model  and  will  be  internally  corrected  if  it 
occurs.   The  data  output  is  still  valid  if  this  error  occurs 
even  though  an  error  message  is  printed  out.   Other  run  time 
errors  or  errors  in  computing  the  radial  turbine  data  should 
be  cause  for  re-examination  of  the  input  data  provided. 
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APPENDIX  J? 
EVALUATION  OF  THE  ZENER  CYCLE 
F.l   General 

The  Zener  Cycle,  as  described,  by  Clarence  Zener  (15), 
is  based  upon  the  therrnocline  existing  in  the  oceans  of  the 
world.   According  to  Zener,  it  should  be  possible  to  take 
advantage  of  this  temperature  difference,  (20  C  is  used  as 
an  example)  to  power  a  thermal  cj^cle  which  could  in  turn 
produce  useful  energy.   This  supposition  is  based  upon  the 
assumption  that  technology  could  produce  the  necessary 
equipment.   Assuming  that  this  is  possible,  Zener  estimates 
that  "the  tropical  oceans  in  the  year  2000  could  supply  the 
whole  world  with  energy  at  a  per  capita  rate  of  consumption 
equal  to  the  U.S.  per  capita  rate  in  1970,  and  suffer  only  a 
one-degree  C  drop  in  temperature".   This  is  estimated  to  be 
about  60  billion  kilowatts.   Zener  proposes  a  thermal  cycle 
using  ammonia  as  the  working  fluid  and  operating  between  5  C 
and  25  C.   His  estimated  plant  size  is  about  8000  cubic  feet. 

The  cycle  consists  of  taking  high  pressure  liquid  at 
25  C,  expanding  it  through  a  turbine,  and  condensing  the 
vapor  at  5  CJ .   The  fluid  would  then  be  pressurized  and 
heated  to  25  C  to  begin  the  cycle  again.   Sea  water  would 
be  used  to  both  heat  and  cool  the  fluid  through  the  use  of 
heat  exchangers.   Ammonia  is  suggested  as  a  working  medium. 
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P.  2   Working  Fluid 

Without  exploring  the  heat  exchanger  problems  involved 

in  detail,  the  feasibility  of  the  turbines  needed  for  such 

a  cycle  was  evaluated.   The  calculations  were  based  upon  a 

plant  output  of  approximately  100,000  kilowatts  with  an 

assumed  turbine  efficiency  of  90  percent.   Based  upon  this 

the  flow  rate  for  a  specific  fluid  was  calculated.   Ammonia 

and  Preon  21  were  chosen  as  the  fluids;  ammonia  because  it 

was  proposed  by  Zener,  and  Preon  21  because  it  is  typical  of 

"the  recently  developed  refrigerating  fluids"  mentioned  by 

Zener.   The  cycle  operating  points  are  summarized  as  follows, 

along  with  the  necessary  mass  flow  rates. 

Point         Enthalpv  Temp  Pressure 

(BTU/lbm)  (°R)  (psia) 

AMMONIA 

1  628.74  528  124.4 

2  612.575  510  89.19 
2S           610.779            510  89.19 

PREON  21 

1  127.56  528  22.2 

2  124.323  510  15.33 
2S            123.963  510  15.33 

MASS  FLOW  KATE 
Ammonia  -  6050  lbm/sec 
Preon  23  -  30220  lbm/sec 
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F.3   Range  of  Input 

The  inputs  chosen  for  the  radial  turbine  were  based 
upon  a  back  calculation  from  specific  speed  to  RFM  and  number 
of  units  in  parallel.   The  K.PM  was  quantified  as  being  3600/N 
where  N  is  an  integer  number.   This  was  done  because  it  is 
assumed  that  the  machine  would  be  used  to  produce  60  hertz 
power  without  benefit  of  a  reduction  gear.   The  number  of 
units  was  then  specified  to  provide  a  variety  of  machines 
covering  the  range  of  specific  speed  recognized  as  feasible, 
section  2.4.9.   Twelve  rotor  blades  were  specified  for  the 
ammonia  trials  and  15  for  Freon  21.   This  brackets  the  usual 
design  range. 

The  axial  turbine  inputs  were  based  on  combinations  of 
parameters  used  to  achieve  high  efficiency.   A  large  slow 
machine  was  favored  with  high  reactions  and  few  stages. 
This  was  not  immediately  apparent  at  the  outset,  but  after 
several  runs  with  no  combinations  yielding  a  feasible  turbine, 
it  was  indicated. 

A  summary  of  design  inputs  for  one  run  is  shown  below. 

The  feasibility  column  indicates  whether  the  design  meets 

the  feasibility  criteria  of  the  model,  sections  2.3.8  and 

2.4.9. 

Axial  Turbine  Inputs-Ammonia 

it    Stages     Reaction       UPM      Pf        Units       Feasible 
,  lm 

2  .5  900       60         4  Yes 


§    Stages     Reaction 


10 

5 

2 
10 

5 

2 
10 

5 

2 

5 
10 

5 
10 

2 

5 
10 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 
20 
20 
10 
20 
10 
20 
20 
20 
20 
20. 


2 
5 
10 
2 
5 


.5 
.5 
.5 
.5 
.5 
.5 
.5 
.5 
.5 
.5 
.5 
.5 
.5 
.5 
.5 
.5 
.5 
.5 
.6 
.7 
.8 
.5 
.6 
.7 
.8 
.5 
.6 
.7 
.8 
.2 

•  c 

.2 

O 

•  <~ 

.3 
.3 
.3 
.3 
.3 
.3 


.7 
.7 
.7 
.7 
.7 


Axial 


RPM 

^m 

720 

60 

720 

60 

720 

60 

900 

60 

900 

60 

900 

60 

720 

60 

720 

60 

720 

60 

900 

'  60 

900 

60 

720 

60 

720 

60 

720 

60 

900 

60 

900 

60 

900 

60 

900 

60 

900 

60 

900 

60 

900 

60 

900 

70 

900 

70 

900 

70 

900 

70 

900 

80 

900 

80 

900 

80 

900 

80 

720 

60 

600 

60 

600 

60 

600 

50 

600 

60 

600 

60 

600 

50 

600 

40 

600 

30 

600 

20 

Turbine 

Inputs 

600 

60 

600 

60 

600 

60 

720 

60 

720 

60 

Units 

4 

4 

4 

2 

2 

2 

2 

2 

2 

4 

4 
10 
10 
10 
10 
10 
10 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 


20 
20 
20 
20 
20 
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feasible 

No 

Yes 

Yes 

No 

No 

Yes 

No 

No 

Yes 

Yes 

No 

Yes 

Yes 

Yes 

Yes 

Yes 

No 

Yes 

Yes 

Yes 

Yes 

No 

No 

Yes 

Yes 

No 

No 

No 

No 

No 

No 

Yes 

No 

Yes 

No 

Yes 

Yes 

Yes 

Yes 

Yes 
No 
No 
No 
No 
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#  Stages  Reaction  RPM     °^  Units       Feasible 

°  lrn 

10  .7  720       60  20            No 

2  .7  720       70  20            No 

5  .7  720       70  20            No 

10  .7  720       70  20            No 

2  .7  720       80  20            No 

5  .7  720       80  20            No 

10  .7  720       80  20            No 

2  .2  720       60  20            No 

2  .3  720     .  60  20            No 

2  .4  720       60  20            No 

2  .6  720       60  20            No 

2  .7  720       60  20            No 

2  .8  720       60  20           Yes 

2  .9  720       60  20            No 

2  .7  600       70  20            No 

5  .7  600       70  20            No 

10  .7  600       70  20            No 

2  .7  600       80  20            No 

5  .7  600       80  20            No 

10  .7  600       80  20            No 

Radial  Turbine  Inputs-Ammonia 

RPM  Units  #  Blades  Feasible 

3600  13  12  Yes 

3600  12  12  Yes 

3600  11  12  Yes 

3600  10  12  Yes 

3600  9  12  Yes 

3600  8  12  Yes 

3600  7  12  Yes 

3600  6  12  Yes 

3600  5  12  Yes 

3600  4  12  No 

3600  14  12  Yes 

3600  15  12  Yes 

3600  16  12  Yes 

3600  '  17  12  Yes 

3600  18  12  Yes 

3600  19  12  Yes 

3600  20  12  Yes 

3600  21  12  Yes 

3600  22  12  Yes 

3600  23  12  Yes 

'3600  24  12  Yes 

1800  1  12  No 


RPM 

1800 
1800 
1800 
1800 
1800 
1200 
1200 
900 


3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 
3600 


Units 

2 
3 

4 
5 
6 
1 
2 
1 


#   Blades 

12 
12 
12 
12 
12 
12 
12 
12 


Radial  Turbine  Inputs-Freon  21 


250 
500 
750 
1000 
1250 
1500 
1750 
2000 
2250 
2500 
2750 
3000 
3250 
3500 
3750 
4000 
4250 
4500 
4750 
5000 
5250 
5500 
5750 
6000 
6250 


15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
15 
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Feasible 

Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 


No 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 


.F.4   Sample  Output 

The  output  for  both  working  fluids  produced  many  designs 
which  could  be  called  feasible.   A  representative  sample  of 
axial  and  radial  designs  for  each  fluid  is  included  on  the 
following  pages.   These  designs  are  not  necessarily  optimums, 
as  the  object  of  the  runs  was  to  test  the  program  rather 
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than  optimize  a  design.   Data  from  these  runs  was  used  to 

substantiate  the  evaluation  of  the  model,  sections  3.5  and 

3.6.   These  examples  do,  however,  illustrate  several  points 

about  the  Zener  Cycle  and  will  be  the  basis  for  the  next 

section . 

Axial  Turbine-Ammonia 

Assumed   5.  Stages   0.50  Reaction  at  mean   720.  RPM 
Assumed   60.00  ALFAMN   4.  Units   6050.00  m 

Last  stage  tip  rel  Mach  no.   0.306  Last  stage  hub  real  hach 

no.   0.314 
First  stage  tip  diam.  8.51   Blade  height  0.9599  it    blades  50. 

Lff.  0.9589 
Last  stage  tip  diam.  8.84   Blade  height  1.2949   #  blades  38. 

Eff.  0.9598 
Machine  efficiency   0.9594 
Max  cent,  stress   3003.50 
Unit  length   2.492 
Aspect  ratio,  first  stage  2.17   Last  stage  2.34 

Velocity  Triangles 

First  Stage 

Hub 

UHN  248.3  VTHN  325.9   VX  164.2   VHN  365.0   VHN  181.7 

WHR  297.7  ALKAHN  63.3   BETAHN  25.3   BETAHR  -56.5 

Mean 

UMN  284.5  VTMN  284.5   VX  164.2   VMN  328.5   VMN  164.2 

WMR  328.5  ALFAMN  60.0   BETAMN  0.0   BET AMR  -60.0 

Tip 

UTN  320.6  VTTN  252.4   VX  164.2   VTN  301.1   WTN  177.9 

WTii  360.3  ALFATN  56.9   BETATN  -22.6   BETATR  -62.9 

Last  btage 

Hub 

UHN  2  3  5.6  VTHN  343.4   VX  164.2   VHN  380.6   WHN  196.4 

WHR  287.2  ALFAHN  64.4   BETAHN  33.3   BETAHR  -55.1 

Mean 

UWN  284.5  VTMN  284.5   VX  164.2   VMN  32  8.5   WMN  164.2 

WMR  328.5  ALFAMN  60.0   BETAMN  0.0   BETAMt  -60.0 

Tip 

UTN  333.3  VTTN  242.8   VX  164.2   VTN  293.1   WTN  187.5 

WTR  371.5  ALFATN  55.9   BETATN  -28.9   BETATR  -63.8 
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Axial  Turbine-Freon  21 

Assumed   2.  Stages   0.70  Reaction  at  mean   720.  RPM 
Assumed   60.00  ALFAMN   20.  Units   30220.00  m 

Last  stage  tip  real  Mach  no.   0.741   Last  stage  hub  real 

Mach  no.  0.542 

First  stage  tip  diam.  8.87   Blade  height  1.9793   #  blades  84. 

Eff.  0.9712 

Last  stage  tip  diam.  9.67   Blade  height  2.7744  #    blades  23. 

Eff.  0.9701 

Machine  efficiency  0.9706 

Max  cent,  stress  5878.27 

Unit  length  0.802 

Aspect  ratio,  first  stage  10.00  Last  stage  4.59 

Velocity  Triangles 

First  Stage 

Hub 

UHN  185.2  YTHN  218.7   VX  90.0   VHN  236.5   WHN  96.0 

WHH  205.9  ALFA1IN  67.6   BETAIIN  20.4   BETAHR  -64.1 

Mean 

UMN  259.8  VTMN  155.9   VX  90.0   VMN  180.0   VMN  137.5 

WMH  275.0  ALFAMN  60.0   BETAMN  -4  9.1   BETAMR  -70.9 

Tip 

UTN  334.4  VTTN  121.1   VX  90.0   VTN  150.9   WTN  231.5 

WTIi  346.3  ALFATN  5  3.4   BETATN  -67.1   BETATR  -74.9 

Last  Stage 

Hub 

UHN  155.2  V'THN  260.9   VX  90.0   VHN  276.0   WHN  138.8 

WHR  179.4  ALFAHN  71.0   BETAHN  49.6   BETAHR  -59.9 

Mean 

UMN  259.8  VTMN  155.9   VX  90.0   VMN  180.0   WMN  137.5 

WMR  275.0  ALFAMN  60.0   BETAMN  -49.1   BET  AMR  -70.9 

Tip 

UTN  364.4  VTTN  111.2   VX  90.0   VTN  143.0   WTN  268.8 

WTIi  375.4  ALFATN  51.0   BETATN  -70.4   BETATR  -76.1 
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Radial  inflow  Turbine-Ammonia 

Assumed   3600  RPM   17.  Units   12.  Blades   6050.00  WW 

Nozzle  inlet  diam.  6.99   Blade  height  0.2338 

Nozzle  outlet  diam.  4.63 

Rotor  inlet  diara.  3.70   Outlet  tip  diam.  2.91 

Rotor  hub  diam.  1.52 

Unit  length  1.09 

Max  scroll  area  3.972 

Specific  speed  0.127727 

Efficiency  0.9398 

Velocity  Triangles 

Nozzle  Inlet 

VI  362.8   VT1  229.8   VM1  280.7   ALFA1  39.3 

Rotor  Inlet 

U3  696.8   VT3  580.7   VM3  364.9   V3  685.8   W3  382.9 

ALEA3  57.9   BETA3  -17.7 

Rotor  Outlet 

Hub 

U4II  287.4   V4H  241.4   W4H  375.3   BETA4H  -50.0 

Mean 

U4M  418.1   V4M  241.4   W4M  482.8   BETA4M  -60.0 

Tip 

U4T  548.8   V4T  241.4   W4T  599.5   BETA4T  -66.3 
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Radial  Inflow  Turbine-Ereon  21 

Assumed   3600  HPM   1250.  Units   15.  Blades   30220.00  WW 

Nozzle  inlet  diarn.  2.75   Blade  height  0.1046 

Nozzle  outlet  diam.  2.04 

Rotor  inlet  diam.  1.62   outlet  tip  diam.  1.23 

Rotor  hub  diam.  0.72 

Unit  length  0.45 

Max  scroll  area  0.628 

Specific  speed  0.113886 

Efficiency  0.9345 

Velocity  Triangles 

Nozzle  Inlet 

VI  173.4   VT1  124.8   VM1  120.4   ALFA1  46.0 

liotor  Inlet 

U3  305.8   VT3  265.0   VM3  128.1   V3  294.3   W3  134.4 

ALEA3  64.2   BETA3  -17.7 

Rotor  Outlet 

Hub 

U4I1  135.1   V4H  105.9   W4H  171.7   BETA4H  -51.9 

Mean 

U4M  183.5   V4M  105.9   W4M  211.8   BETA4M  -60.0 

Tip 

U4T  231.8   V4T  105.9   W4T  254.9   BETA4T  -65.4 
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F.5   Comment  on  the  Zener  Cycle 

Using  the  sample  output  shown  in  the  previous  sections, 
it  is  possible  to  calculate  a  volume  required  for  each  machine 
Multiplying  this  by  the  number  of  units  gives  a  volume  to 
produce  100,000  kilowatts  of  power.   This  approximate  volume 
is  600  ft3  for  Axial  with  Ammonia,  1200  ft3  for  Axial  with 
Preon  21,  700  ft3  for  Radial  with  Ammonia,  3700  ft3  for 
Radial  with  Freon  21.   This  is  interesting  in  that  it  shows 
ammonia  to  be  the  more  advantageous  of  the  two  fluids  in 
terms  of  size.   This  would  be  expected,  however,  in  that  for 
the  cycle  points  defined,  ammonia  has  a  greater  enthalpy 
change  per  pound  than  Freon.   i'reon  has  the  advantage  of 
being  safer  to  work  with,  however. 

A  second  fact  is  apparent  in  the  volume  consideration 
in  that  the  estimated  volume  was  in  the  order  of  one  to  two 
thousand  cubic  feet.   Zener  estimated  a  turbine  module  to  be 
8  x  8  x  40  feet  or  about  2560  cubic  feet.   This  shows  that 
Zener  '  s  estimation  was  reasonable  for  plants  of  100,000 
kilowatt  output.   Assuming  this  size  plant,  then  to  supply 
the  output  of  60  billion  kilowatts  mentioned  in  section  F.l, 
600,000  individual  plants  would  be  required.   This  also 
assumes  that  the  heat  exchangers  are  of  comparable  size. 

Comparisons  could  be  made  for  the  amount  of  working 
medium  required,  amount  of  cooling  water  and  others.   That 
is  not  really  the  point  of  the  evaluation,  however.   Lxami- 
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nation  of  the  sample  designs  does  show,  for  at  least  the 
turbines,  that  the  proposal  is  at  least  within  the  realm  of 
possibility  if  not  reasonability . 
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